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XXIV. BAND FUNFTES HEFT 1956 


Formeln fiir die Massenkriafte und kinematischen Zusammenhinge 
bei geschrankten Schubkurbelgetrieben 


Von E. Mewes 


1. Allgemeine Beziehungen. Fiir das in Abb.1 skizzierte System wird der Verlauf der Koor- 
dinaten x, y und ¥ mit dem sich stindig andernden Kurbelwinkel y verfolgt. Enthalten sind 
dabei auch folgende Sonderfalle: Fiir € =1 und 7, =0 (Sonderfall 1) ergeben sich die Wirkungen 
hin- und hergehender Massen. Fiir  =7 =0 (Sonderfall 2) ergeben sich die Wirkungen unwuchtig 
rotierender Massen. 

Wenn man von dem bei ungleichférmigen Umlauf der Kurbel um die Kurbelwellenachse 
drehenden Moment der Massenkrafte der umlaufenden Massen absieht, so kann man bei un- 
wuchtig umlaufenden Massen diese auf den AnschluBpunkt 1 (Abb.1) reduzieren. Andererseits 
kann fiir die Ermittlung der Massenwirkungen die Masse des hin- und hergehenden Teils im 
Anschlu8punkt 2 angesetzt werden. Mit diesen Massen und der Masse des Verbindungsgliedes 
(Koppel) zusammen entstehe der Massenmittelpunkt (Schwerpunkt) S, der sich mit der Koppel 
starr verbunden bewegt. 


Abb. 1. Geschranktes Schubkurbelgetriebe mit einseitiger Koppel (Pleuelstange). 


Der Fall 740 kommt z.B. bei V-Motoren an Hauptpleuelstangen vor, an denen jeweils 
ein Nebenpleuel angelenkt wird. Auch bei den sog. Horden- und Ladenschiittlern in Dresch- 
maschinen und anderen Siebvorrichtungen kommen Systeme gemafs Abb.1 vor, wenn man bei 
der Bahn des zweiten AnschluBpunktes der Koppel von den geringen Abweichungen von der 
Geraden absieht, die bei Halterung durch lange Lenker auftreten. 


Allgemein ist nach Abb.1: 
x =reosy + élcosy+ylsiny, | 
y =rsiny—élsiny + yl cosy, J 
also mit 
A= Til, 


é ‘jee 
ee or xs | 


~|8 


(2) 


oe 


—— ae eT | 
:. SS) — er sin ¥ 7 cosy. 


A 
Zum Unterschied von Zwischengleichungen sind bei den Gleichungen, die fiir die Ergebnisse der 
Schwerpunktkoordinaten und Massenkrafte und -momente benétigt werden, die Gleichungs- 
nummern unterstrichen worden. 

21 
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Im Sonderfall 1 ist 
By 1 
te aay ae I 7, | | 


ae cle =| 
Fah iak torte eke) Soman ae 


mit B =b/I. Allgemein gilt aus r sin y —I sin y = b (nach Abb.1, vgl. auch Schrifttum!) | 
siny =Asiny — 6. (4) 
Dabei ist sin y in einfacher Weise von y abhangig: anders ist es fiir cos y und y. | 
Aus (4) erhalt man durch Ableitung nach der Zeit mit w =y 
x cosy =Aw cosy, 


also | 
y = ho coy" (5), 

Weiter ist | 
cos y = 1 — sin? y = [1 — (Asin y — f)?]!”. (6) | 


Da in y auch cos y mit dem Exponenten —1 vorkommt, wird allgemein untersucht fiir be- 
liebiges m 
cos™ ¥ = [1 — (Asin y — f)?]™”. (7). 
Auf dem Wege zu Entwicklungen in Fourier-Reihen wird dies zundchst nach dem binomischen | 
Satz in eine (unendliche) Potenzreihe von ( sin y entwickelt. Aus } 


cos™ y == (l— f? + 2B Asiny — 2? sin? p)ymi2 


m/2 
= (1 — fp)? (1 tay os Asin y — iF 2? sin? ¥) (8) 
wird erhalten 
Asi Asin py \2 
cos” y= (1— BA) fay + a Te + my (Tem) + -- | (9) 
mit 
ea 1, 
2 
a, = ei 2 i) 5 
2 2 
a= ("5 ate jap, 
: : (10) 
m= ("))8—("5)(1)260—P). 
m/2 2 2\ /2 
a =(")) 166 —("P)() 4 06) +("P)(3) a—pr. 
(Fiir 6 = 0 vereinfachen sich die Beziehungen erheblich.) 
Setzt man statt (9) 
cos™ y = by + bjAsiny + b,/? sin?y +--- (11) | 
dann gilt (nachweisbar nach dem Satz von Taylor, demzufolge 
oa : ey os of(—B) Asin p 0 f(— B) A? sin? p | 
S(—B +Asin y) =f ( B) 4 O38). malt : a(— B)? 2! ream (12) 
SE) titer en 
ped ih Oa | 
bn ae apy’ (13) 
also auch (als Rekursionsformel) | 
ip eee Se ae Obn—1 | 
Maa n Op (14) 
* Z.B. C, B. Biezeno u. R. Grammel, Technische Dynamik, Bd. Il, Kap. XI, 2. Aufl. Berlin 1953. 
| 
t 
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Es ist 
by = (1— p22, 
b, = (1 = pe) ay, (15) 
Nach (9) und (11) besteht die Beziehung 
b, = —— (1 =e azar gy Qn » 16 
-also auch P) oy) 
b, a (1 f= , pA)sia— 4-1 Goi (17) 
Aus (14) wird mit (17) entwickelt 
L [f/m = aa, 
b, = — n (5 n as 1) (1 — p?)™ : (— 2 B) an—) + (1 — B?) Jizan) ap eh) 


also mit (16) 


1 
a, = — 
n 


26 & pal 1) aa Bae 3) oe | (18) 
Mit dieser Rekursionsformel kann stets aus dem vorhergehenden Koeffizienten a,_; der folgende 
Koeffizient a, berechnet werden. 

Es ist z. B. nach a, = 1 
a, = mp (19) 


und 


= = m [(m — 2) — (1—P*)], (20) 


was dasselbe ist wie in (10). Bei Zusammenfassung der Anteile mit gleichen Potenzen von B 
wird daraus erhalten 


1 
ee Pt (1 am) 6]. (21) 
Wir untersuchen weiterhin die Zusammenhange fiir einzelne m. 


2. Unterlagen fiir Lagekoordinaten und Massenkraftkomponenten. Fiir m=I1 erhalt man 
entsprechend (9) 


A Asi 2 
cos Y = (1 — fp)? B ee ; = [ere ( = fs | (22) 
mit ; 

€ dl 
Cn B . 

it 
C5 —, —— 5 2 

1 

Co al ’ 

1 1 
Ds Er mee 


mit der Rekursionsformel 


Gin) Bec Sh) | (23) 


aa 


On 


Setzt man an 
cosy = fy + fi sin y — f, sin? y + f, sin? y — f, sinty + —-°°, (24) 


Sm Ture! Hie Ge Iw, Bip 
1—2n 


fo =(1—B?) ® -e,d"(— 1)" (25) 
21* 
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und mit (23) oder den aus my zu entwickelnden Beziehungen 


und 


en = 28 ae) (a) eee Oe at ae en. 


erhalten 
h=0—pye, 
fi =(L—P)~'? Ba, 
fi = qd — py 2 
fp= (1 — BP BR, 


RRS a: mit g,=1+4f%, 


fe = (1— PI) SB Bs ae 

ee =1+4 126 + 8ft, a 
fy = (1B? BT er ee 

fy = — BY 55 Bs ge = 1424p + 48 ph 4+ SB, 

f= —P) BP gs = 18 + eH + BY, 

fos (lh) tee ce Bio = 1 + 40 B® + 160 pt + 128 6° + = ge, 


Es sei besonders darauf hingewiesen, daB die Koeffizienten bei den Potenzreihen geschlossene 
Ausdriicke sind; jeder der bei (28) aufgeschriebenen Ausdriicke ist vollstandig. 


Fiir Berechnung der Koeffizienten f kann auch verwendet werden: 


(Wee oi! ey) Be + Core: ile : p? ai “een 
apy =1— (71) + Ge on 
ae Gaim ee 


Aus (24) wird entwickelt mit 


ap 1 
sintey = e(),) + gies D1 UM," ) e082 m 9 (30 


a 1 a 9 (PDH IIN 
sin- : Ve g2n—2 > (== 1)" eee) ah (2 Uk 1) a (31) 
m=1 


die Fourter-Reihe 
cosy = E, + F, siny + E, cos 2y — EB; sin3 yp — E, cos4y + +——::: (32) 
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mit den unendlichen Reihen 
i 


Ey = fo ae 
B= f ala Brae Wet 4% 


= i 1 
i= da epee? Sembee 
E : es = 
= my eh Pale 
= 1 
Ey “6 oe Si ae ? 
Fir n=], 2, 3, ... gilt also allgemein 
<a 2k 
E, = 2 wri ne) fH e (34) 
Fir Sonderfall 1 ist nach (3) 
x af : 
— = cosy + = (Eo + E, siny + E, cos 2 yp — E, sin 3 y — E, cos4y 4+ + ——-::). (35) 
Man schreibt dafiir auch 
x : 1 il ; 
— = Ay + cosy + Ay sin y + Ase0s 2 —> As sin 3p —= A, cos 4 + +——- = (36) 
oder 
4 + cosy + B, si ae Z d i } 
A, y 1 Siny +> B, cos po, By sin 3 yp Be cos Aye jt (37) 
mit Ay = Bia und fiir i = Ih, By Shoo 
A=) Doss Led. (39) 
Fur den allgemeinen Fall nach Abb. 1 gilt nach (2) mit (4) und (32) 
= = cosy +n (siny 1) + 7 (Bo +B; sin y + By cos 2p — Fy sin 3 y 
/ —— Ecos 44 -- + +), (40) 
= siny —E(siny— 7) + [Bo + Bi siny + E, cos 2y— E, os 
— E,cos4y ++——°-:').J 
Daraus ergibt sich 
~ = —siny + (7 + €B,) cosy—E CEST Sapa i 
(41) 
7 =(1—£ + 7B) cos yp — n B, sin 2 py — 7, Bz cos 3p + +——""' | 
und fiir m = konst. 
—— =cosy + (1 + &A,)siny +A, cos 2y—EAysin3y —EA, cos dy + +——""s| 
= (42) 
—-, =(1—é+y A,) sin y + 7A, cos 2y—n A, sin 3y —7 A, cos4ap a= = | 


HeiBen die Komponenten der freien Massenkraft X und Y, so sind die in (42) dargestellten 
GréBen auch X/mrw? und Y/mrw?, wenn m die (in S zusammengefabte) bewegte Masse ist. 


3. Drehungen des Verbindungsgliedes (Koppel, Pleuel) im Schubkurbelgetriebe. Der Verlauf 


von x ist noch zu bestimmen. Bei dem in (5) vorkommenden I/cos y ist gemaB (9) fir m= — 1 


iI : . Asin y . {Asin 
Sd pose, 4-4, ee sent (TF) eee | (43) 


cos ¥ 
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mit ip = 1 Tiel shbie 7p == IL; M4, Bb, ood | 
a < a 1/2 k c—n p n—k 
UDI. = > (— LE ( k ( Te x (4 py (1 — By? v (44) 


tan—1 = 28 | 3 (Gar ieee ee (4p ee (45) 


k=n+1 


bzw. mit der Rekursionsformel nach (18) 


fee p (2p eee ae | (46) | 


op 


Die Koeffizientenanteile sind zweckmaBig jeweils nach Potenzen von fp zu ordnen. Damit wird | 
z. B. 
i = is 
: 1 
eae + 6, 
7 3 
i= —B(z +B). 


Bei Ansatz von 


ay =k, —k, siny + ky sin? y — kz sin? yp + —:-- (47) | 
ist 
ky = (1 — B?)—? (48) 
und 
key = (— 1)" (1 — fy" "3,2", (49) | 
also “a 


ky == (L— 6)" B27, 
ky =(1—B)-39 5 1, mit _ = 14262, 


(50) 
x @ : 2 
ky = (1p) 36 mit = 14562 
usw. zu berechnen aus (49) mit (46). Fiir cos y sin" y erhalt man infolge der Beziehungen (30), 
(31) und 
: Laer ; 
cosy sin my = > [sin (m + 1) y + sin (m— 1) y], (51) 
l 
cos py cos my = — [cos (m + 1) y + cos (m— 1) y] (52) | 


die Beziehungen 


cos p sin?” y = ae [> (— 1)**! ey, cos (2 k — 1) y + (—1)" cos (2 n + 1) v| : (53)5 | 
a | 
cos py sin?™—] y = as bs (— 1)*+!d,, sin2ky + (—1)"—-'sin2n v| (54) | 
= a | 
a ‘ ae 4. = boa as <P ee .) tO eae o eo) 
und 
d,, = Eee i (. ea) a — if eC fir n—k>0. (56) 


Damit ergibt sich die Fourier-Reihe (vgl. Biezeno-Grammel) 


cos y/cos y = D, cos yp — D, sin 2 y — D, cos 3 p + +,——---- (57) 
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| 


mit den unendlichen Reihen 


1 1 
D, = hy oh Ky + 7 Cas Mis ee, 


it I 
D; = a Beier || ae) 


1 1 il 
D, Sg ar oy ks gg dey he es 


1 1 5 
D, = 8 ks ee | (59) 
Dabei ist z. B. 
He 2 é 4 4 ¢ 6 
-(i)-@)-% #-(A()-# aC) 
4 4 6 6 
a-(i}-(f)- Ge (es 
0), ie 
alice 
3 3 3) 5 
Berm eee a) ig ® 
BY 5 
2-()-()- 
Fir die Drehbeschleunigung 7 gilt (vgl. Biezeno-Grammel) 
— y/Aw? = C, siny + C, cos 2p — C, sin 3y — C, cos4y + +——-::> 
— w/w* (D, cos y — D, sin 2 y — D, cos 3 y + D, sin 4y + —-:->) (60) 
mit ne 
(C5, aif Deg. (61) 
Durch Integration von (5) mit (57) wird erhalten 
y=A(—F, +f siny + F, cos 2 yp — F; sin 3 y — F, cos 4y + +——-::-) (62) 
mit 
7 D,, C, =e 
Dike sate (hie ge I, Py By 6 ove (63) 
und 
il oo 
F, = 7 are sin Bf + > (— 1) PtP (64) 
n=1 


4. Berechnung der Massenkraftkomponenten fiir ein Beispiel. Um die Anwendung einiger 
Formeln zu zeigen, wird ein Beispiel behandelt. Dabei sei die Kurbelwelle dreifach mit je 120° 
Kurbelversetzungswinkel gekrépft. An jeder Kurbel sei ein gleich groBer Schiittler als Koppel 
angelenkt. Die Bahnen der Gleitstiicke liegen nebeneinander. Berechnet werden im folgenden 
die Massenkraftkomponenten der einzelnen ,,Schiittlerladen“ und die freien Massenkraftkompo- 


nenten des Gesamtaggregats. Letztere sind 
X,=3mrw? & (— Ag sin 3 p + Ag cos 6 yp — Ay sin 9 ae (65) 
Y, = 3mro? ny (— A; sin 3 py + Ag cos 6 p — Ay sin 9 y+ —*""). 

Auf die Momente wird dabei hier nicht eingegangen. 

Gegeben ist B = 0,4. Dafiir wird erhalten 
(Lh *) AAO ONT (= 6?) 87 = 1530, el — fp) 9 84, 
(Veja 0 091, (1 — B?)—*? = 1,55, (1 f2)- "2 = 2,19 usw. 
und auBer ¢ = 1, e, =— 0,5, in = 1, e, = 0,4, €3 = 0,2, €, = — 0,205... 22400, = 0,4, 1, == 0,60, 
iz = — 0,664,... fy = 0,917, ky = 1,091. 
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Gegeben ist ferner A = 0,2. Damit wird erhalten 
f; = 0,0873, f, = 0,0260, f, = 0,00248, 7, = 0,000004, J, —= 0,000128 >. 
ky = 0,104 ky = O04. kh, 00098 ee, 


E, = 0,917 — 0,013 — 0,904, 
E,= _ 0,0873 + 0,0019 + 0,0001 — 0,0893 , 
i 0,0130 + 0,0003 = 0,0133 , 
ee 0,00062 + 0,00004 — 0,00066 , 
Ez 0,0001 , 


Ay = 4,52, A, = 0447, Al = 0,260, A, = 0,030, 94, = 0.0016 2 
Bre 0,447, By. =0,133, 8, = 0,010, 08 =) 0004 ee 
Dy = 1,001 =. 0,010 1,101 = «=, os 
D, ~ 0,0520 + 0,0025 = 0,0545, C,—0,109, F, = 0,0273. 


Somit ist x = 0,2 (— 2,08 + 1,10 sin y + 0,03 cos 2y— +--+), 
= age = 1,10 sin w + 0,11 cos 2 y — + 


Gegeben ist ferner als Umdrehungszahl 600 U/min. Damit wird 
Gi 62-8isel Aw — 19 0ise2s 
y = — 870 sin p — 87 cos2p + —---. 
Fiir ein Massentragheitsmoment einer Koppel um die Querachse durch ihren Schwerpunkt von | 
0,06 kgm s® ergibt sich als Moment in kgm von der Massenwirkung einer Koppel um diese Achse | 
Mx 52 siny +5 cos2p—-+-°-- 
Gegeben ist ferner € = 0,4, 7 = 0,1. Damit wird fiir cine Koppel 


= cosy + 0,279 sin y + 0,106 cos 2 y — 0,012 sin 3y—4+::-, 


mr ow 
os a 0,645 sin Y at 0,027 cos 2 ie 0,003 sin 3 yy — + im 
Fur das Gesamtaggregat mit drei Koppeln (Schiittlern) ist 
xX ; Yy ; 
rok = 07036 sin 3 y lpheap oes Sie = — 0,009 sin3y + —-:::. 


Die Koeffizienten bei den weiteren Gliedern sind viel kleiner als die aufgeschriebenen. Fiir 
mr @ = 1000 kg ist X, max = 36 kg, Y,max = 9 kg und der GréBtwert der resultierenden freien 
Massenkraft Rg max = /36? + 9? ~ 37 kg. 

Fiir einen Schiittlerkasten ergibt sich nach obigem als fiir dieses Beispiel (6,4) brauchbare 
Naherung fiir die Massenkraftkomponenten 


X ~ mr wo € cos p 4 rn ea 


ee oe ee 
cial tirdap + Jo 
+r tak aa . Pao a Soe aah 


9 ; 
a ve + a dn Sel 
Solche Formeln, bei denen die Abbruchstellen fiir die Reihen bewuBt an bestimmte Stellen 
gesetzt wurden, kénnen aber fiir andere Beispiele wegen gréBerer Fehler unzweckmaBig sein. Man _ 


sollte dann stets wieder von den vorher angegebenen, vollstindigen Reihenangaben (42) mit 
(39), (33), (25) und (23) ausgehen. 


L -+- os Mayo ++: |! sin wy 


++ 


(66) 


Ya mrot(l1—é + 


Aus dem Institut fiir Landtechnische Grundlagenforschung der Forschungsanstalt fiir Landwirtschaft, 
Braunschweig-Vélkenrode. 
(Eingegangen am 17. Januar 1956.) 
Anschrift des Verfassers: Dr.-Ing. Ernst Mewes, (20b) Braunschweig, Bundesallee 50. 
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Strémung in Blasleitungen 
Von H. Hesse 


1, Einleitung. Bei der pneumatischen Foérderung wird das feste Gut in eine Rohrleitung 
gegeben, in der es von einem Luftstrom mitgerissen und dadurch weiterbewegt wird. Die auf- 
gegebenen Kérner haben zunichst im wesentlichen die Geschwindigkeit null und werden durch 
die strémende Luft beschleunigt. AuBer der dazu nétigen Arbeit mu® von dem Luftstrom noch 
diejenige aufgebracht werden, die zur Uberwindung der Reibung an der Rohrwand dient. Die 
folgenden Uberlegungen befassen sich mit diesen Vorgingen; dagegen soll von denen, die durch 
Steigung oder Fall der Leitung bedingt sind, abgesehen werden!. 


2. Ableitung der Gleichungen. Es seien bezeichnet mit L das in der Zeiteinheit durch die 
Leitung strémende Luftgewicht [kg/s], G das in der Zeiteinheit durch die Leitung geférderte 
Gutgewicht [kg/s], y; das spezifische Gewicht der Luft [kg/m?], y, das spezifische Gewicht des 
festen Gutes [kg/m?], v die mittlere Luftgeschwindigkeit iiber den Rohrquerschnitt [m/s], ¢ die 
mittlere Geschwindigkeit des festen Gutes [m/s], F der Rohrquerschnitt [m?]. 

Fiir die Strémung im Rohr gilt die Kontinuitatsgleichung 


L G 


—F=F,+F,. (1) 


Wenn man sich das Gemisch in zwei Stréme getrennt denkt, stellt der erste Summand den 
Querschnitt F, fiir die Luft, der zweite den fiir das feste Gut F, dar. Es werde ein unendlich 
kleines Stiick des Férderweges s [m] abgegrenzt. In dem entstehenden Raumelement (F ds) 
befindet sich an Luftgewicht F; y; ds und an Kérnergewicht Fy, ds. Nach Gl. (1) ergibt sich 
das Verhaltnis des in dem Raumelement befindlichen Luftgewichtes zu dem des festen Gutes 
wie L zu G v/c. 

Nach dem Impulssatz mu®B der Unterschied zwischen den BewegungsgréBen der in das 
Raumelement in der Zeiteinheit ein- und austretenden Stoffmengen gleich dem Unterschied der 
auf den Ein- und Austrittsquerschnitt wahrend dieser Zeit wirkenden Kraft sein, also gleich 
‘dem Unterschied der Produkte aus dem Querschnitt Ff und dem Druck P [kg/m?]. Durch die 
BewegungsgroBen der Luft und des festen Gutes ist aber der Impuls, der durch die Reibung auf 
das Rohr ausgeiibt wird, nicht erfaft. Er soll daher vor Aufstellung der Impulsgleichung von dP 
abgezogen werden. Man teilt die Reibung in die durch die Luft verursachte dP, und die durch das 
feste Gut verursachte dP,. Es ist 


apa 


v? V1 2 9 
Dg (2) 
worin A die Widerstandsziffer der Rohrleitung gegen Luft und D [m] den Rohrdurchmesser be- 
deuten. Fiir dP, ]a8t sich wahrend der Beschleunigung des Férdergutes die von Barth? angege- 
bene Gleichung verwenden: 


Ge G 
es Bek 3 
dP, = kata (Fe) o- (3) 
G/(F c) ist ein gedachtes spezifisches Gewicht, das entstehen wiirde, wenn man das feste Gut 
gleichmaSig den Raum ausfiillend, iiber den Rohrquerschnitt verteilte. k ist eine Konstante, 
die Widerstandsziffer des Rohres gegen das feste Gut. 


1 Wahrend der Niederschrift dieser Arbeit erschien tiber das gleiche Thema ein Aufsatz von Weidner, 
,,Grundsatzliche Untersuchung iiber den pneumatischen Férdervorgang, insbesondere iiber die Verhaltnisse 
bei Beschleunigung und Umlenkung“‘, Forschung Ing.-Wes., Diisseldorf 1955, S. 145. Da aber in dem vor- 
liegenden Aufsatz immerhin einiges Neue dazu beigetragen wird, besonders was den Vorgang bei groBer Gut- 
beladung und die Aufgabe des festen Gutes ins Rohr betrifft, so diirfte sich seine Veroffentlichung doch noch 
lohnen. 

2 W. Barth, Fortschritte der Verfahrenstechnik S. 63, Weinheim a. d. Bergstr., 1952/53; Weidner, 
Untersuchung iiber den pneumatischen Férdergang, ..., Forschung auf dem Gebiete des Ingenieurwesens 
1955, S. 147. Der dort verwendete Beiwert a, = k/8. 
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Nun lat sich die Impulsgleichung niederschreiben: 
L G 
F (dP — dP, —dP,) = ie dv + = dc 


oder 


L G ds 
UN See umecaren ut 82 Ss iy" (4) 


Diese Gleichung sagt noch nichts aus tiber die Kraft, die zwischen den Bestandteilen L und G 
des Gemisches auftritt. Je nach der Art der Kérner ist diese Kraft proportional der ersten bis 
zweiten Potenz der Relativgeschwindigkeit v — c. Die Abweichung von der zweiten Potenz werde 
mit n bezeichnet. Das kennzeichnende Ma® fiir diese Kraft, die dem Widerstand der Korner 
gleich ist, ist ihre Sinkgeschwindigkeit wy [m/s] in Luft des Betriebszustandes, bei der sich 
Widerstand und Gewicht aufheben. Nach Feststellungen von Barth! andern sich die Sink- 
geschwindigkeit mit der Gutskonzentration G v/(L c) und die Berichtigung n des Exponenten 2 
mit der Anstrémgeschwindigkeit v —c. Das in ¢ Sekunden aufgegebene Gewicht @ t leistet also 
einen Widerstand 


v—c\2—n 
W = 64 re (9) 


Aus einer Impulsgleichung von der Art (4) erhalt man eine fiir die Krafte, wenn man sie durch 
das Zeitelement dt dividiert. Die Kraftegleichung fiir das feste Gut muB noch ein Glied enthalten, 
das den Widerstand des in 1 s aufgegebenen Gutes ausdriickt, weil der Impuls nach (4) auch jede 
Sekunde neu erteilt werden mul}, Sie lautet 
CP ee OP ee 
Maa, = decor gan ar 


oder, wenn man mit Hilfe der Beziehung ¢ = ds/dt das Zeitelement und mit Hilfe der Stetigkeits- 
gleichung (1) F, ersetzt, 


G dP k v—c\2—n G dec 
é 2 ee pee 

Y_ ds pane ! ( ae = ep tds (6) 

Will man die entsprechende Kraftegleichung fiir die Luft aufstellen, so mu man beachten, dab 

dt = ds/c festgesetzt ist. Man erhalt dann 

L dP v yy O(a C\ ee eae, 

cme ta oe | ate rdee (7) 

Wenn man die Kraftegleichungen (6) und (7) mit dem Wege ds multipliziert, so stellen die ein- 


zelnen Summanden abgegebene oder verbrauchte Leistungen dar, Arbeiten in derselben Zeit dt 
sowohl in (6), als auch in (7), da G und L in der gleichen Zeit durch jeden Querschnitt der Leitung 


str6men. Man erkennt, daB die Luft an Leistung ee ds abgibt, aber das feste Gut nur ye ds 


aufnehmen kann. Der Unterschied kann nur in Reibungswarme verwandelt werden. Er ist 


aR =7(2—1) ds = oO eS ie (8) 
t\c c wr” 

3. Integrierung der Gleichungen. Es sei vorausgeschickt, da8 in allen folgenden Ausfiihrungen 
das spezifische Gewicht der Luft y; als unveranderlich angenommen ist. Dann lassen sich unter 
gewissen Bedingungen einzelne Gleichungen oder Teile davon integrieren. 

Erster Fall: F, ist klein gegen Fj. Dann kann man L/(y; v) ~ F setzen. Fiihrt man anstelle 
der Veranderlichen v im Reibungsglied der Impulsgleichung (4) L/(y; F) ein, so 1aBt sich der 
gréBte Teil davon unter Voraussetzung eines gleichbleibenden F geschlossen integrieren: 


ees C AL? kG if 
Ae or ts) stars al Corea) | Jen F2D So) 4 2 DF | cds. (9) 


So 
Die GréBen mit Index 0 bezeichnen den Anfangswert der Integrierung. Zur vollstandigen Lésung 
brauchen wir eine Beziehung zwischen c und s. Dazu benutzen wir am besten die Kraftegleichung 
(6), in der wir das erste Glied vernachlassigen kénnen, wenn F’, vernachlassigt werden kann. 


W. Barth, Chemie-Ingenieur-Technik 26 (1954), S. 31. 
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Dann bleibt v angenahert konstant. Es ergibt sich der Differentialquotient 


ds ~, I 
de. & crs 5 : (10) 


c 


Diese Gleichung ]aBt sich graphisch integrieren, indem man ds/de iiber c auftragt und die ent- 
standene Flache planimetriert. Dann zeichnet man ein Schaubild ¢ iiber s und erhilt wieder durch 
Planimetrieren |¢ ds in (9). 


Fiir n = 0 1aBt sich (10) auch geschlossen integrieren. In die Gleichung 
cdc 


aed hve 
wy wi Je: (‘ STB 


fiihrt man ein 
gvjwp=a und g/we —k/(2D) =h. 
Dann ergibt sich 


av—2ac+hc? ce a/2e is a(c + cy —v) —(c— cq) 0,5 avk/D—heyc 


1 
S— Sy ~ In ———— 
2h av—2Zacq+the a (c +c —v) + (ec —e) 0,5 avk/D —heyc 


| (11) 


Unter der Bedingung n = 0 1aBt sich Gl. (8) teilweise geschlossen integrieren, wenn man sie in 
einzelne Summanden auflést: 


C s s s 
R— R, wal |e 3v2 (s—so) +30 [eds— fe “9 
ei Cc ‘ 
Der Reibungsverlust innerhalb des Gemisches ergibt sich daraus zu 


eR ah een ee 
w h 2avk a(c + ¢)—v) + (ec —e) V0,5avk/D—hc cy 


(12) 


2h AOS eer Wlaee 


So So 


3 2 = 2 : ‘ j [ 
sy Feat ot at sce ee | 


Die beiden letzten Glieder mu man wieder graphisch integrieren. In Abb. 1 sind die Werte von 
c und P uber dem nach der Formel (11) berechneten Weg s fiir Weizen dargestellt. 

Zweiter Fall: Beriicksichtigung von F,. Es gibt auch Férderleitungen, bei denen man F’, 
nicht vernachlassigen darf, z.B. solche fiir Kohlenstaub. Da kommt es vor, dai das Staub- 
volumen (G/y,) gréBer als das Luftvolumen (L/y;) ist®. Ist v konstant, so muB der Rohrquer- 
schnitt F nach der Stetigkeitsbeziehung in der Stromrichtung abnehmen. Eine Losung kann 
fir konstantes v angegeben werden, wenn es angangig ist, fiir F und D in den Gliedern fiir die 
Rohrreibung der Gl. (4) Mittelwerte F,,, und D,, einzusetzen. Das Glied fiir die Beschleunigung 
des festen Gutes formt man mit Hilfe der Stetigkeitsgleichung (1) um. Fiihrt man fiir die Gut- 
beladung G/L die Bezeichnung y ein, so ergibt sich 


c s 
AP = i ee pA ST (s $) Preis) oie 
J ele * 7) : 
Setzt man im ersten Integral als neue Veranderliche 
c il 
a yy v Vg 


1 Siehe Hiitte I 1949, S. 107, Nr. 9und 10. Von denim Logarithmus des zweiten Summanden angegebenen 
Vorzeichen kommt hier nur das positive in Frage; fiir cy = 0 ist der damit versehene Bruch positiv, wenn 
v—c\? k 
We 2gD 
2 Erk, Chemieingenieur 1 (1933) S. 140. 


c2. Das ist immer der Fall. 
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cin, so lat es sich losen, und man erhalt 


s 
1 P Soll hoy tcovg|, vy arose 2kG ARE Qa 
eis Hg AC) scale yee vy + CYg + Ase, | 0+ Dix Ca 
So 
ic p 
73 
— | TL eccanstezesemisek | eee oa ee 
= 
20 S 
[rm/sek] = 


15 


10 


0 & 10 6 20 25 {m] 3% 


Abb. 1. Zunahme der Gutgeschwindigkeit c und Druckverlauf Pin der Rohrstrecke s fiir Weizen (n =0) bei folgenden Annahmen: 
v = 30 m/s = konst., 10 = 9,8 m/s, D= 0,113 m, p = 4, k = 0,0027, y_ = 1260 kg/m$, wy 1,58 kg/m? ~konst., 4 = 0,025, 
P, = 13000 kg/m? abs., c,= 0. 


x7000\g/mt abs > & 


~ 
Ss 
a 


Wegs 


0 5 


70 15 20 


ob 30 (em] 35 
Abb. 2. Zunahme der Gutgeschwindigkeit c, Druckverlauf P und Rohrdurchmesser D fiir einen Staub (Kohlenstaub) bei folgenden An- 
nahmen: Korndurchmesser 0,075 —0,08 mm, wr = 0,217 m/s (Z. VDI 1928, S. 947), x = 3000, v = 5 m/s = konst., Yo = 1400 kg/m’, 
y= 1,75 kg/m* ~konst., k = 0,01, 


= 0,025, n = 0,6 im Durchschnitt gesetzt, zunehmend mit c von etwa 0,2 bis 0,85, berechnet mit 
Hilfe der Oseenschen Gleichung, P, = 15000 kg/m? abs., c, = 0. 
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Weil man in der Kraftegleichung (6) den Summanden Ss iS nicht mehr weglassen darf, ver- 
Ss 


einigt man (6) und (7), um dP zu entfernen und erhalt 
OS iee 1 
dc & ( YI Uv [ = ee De n A 7 » ie eee a! (10a) 


1 pe ie eS 
c c +r wr 8Ly,c 2D, / 


Vv 
5 


Je gréBer die Gutbeladung, um so kleiner ist die durch die Luft verursachte Rohrreibung gegen- 
tiber der durch das feste Gut hervorgerufenen. Bei groBem « wird man meist 2 = 0 setzen 
kénnen. In Abb. 2 ist c, P und D iiber s fiir Staub bei groBer Gutbeladung yw dargestellt. Es 
ist dabei an Kohlenstaub gedacht. Die Beziehung von c zu s wurde nach (10a) durch graphische 


Integrierung gefunden. Auf diese Weise kann man auch die Reibung innerhalb des Gemisches 
nach (8) bestimmen. 


4. Der Beharrungszustand » — konst., c = konst. Im Falle daB F, sehr klein ist gegen F, 
_ kann man in (6) den ersten Summanden wieder gleich null setzen, und es wird 


Echs s ie _- Aes 

2¢D \ Ue ; 
weil die Beschleunigung des festen Gutes auf der rechten Seite der Gleichung verschwindet. Ist 
n = 0, so ergibt sich daraus das von Barth! angegebene VerhAltnis 


Cay 1 
VL + Vk[(2g D) 


(13) 


Darf man n nicht auBer Acht lassen, so lautet die Gleichung 
1 


Veet Sy 


Diese wird man anwenden, wenn n naher an null als an | liegt. Liegt dagegen n naher an 1 als an 
null, so rechnet man vorteilhafter nach der Gleichung 


lee 2gD Ce | ee C/O) gD ise ies 
v ~ Ve wi” (1 7 kv" wa” ko" we" ies (13b) 


Man kann aber auch (1 — c/v)’—" in eine binomische Reihe zerlegen und erhalt dann fiir c/v eine 
Gleichung zweiten, dritten oder héheren Grades je nachdem, bis zu welcher Potenz von ¢/v man 
die Glieder der Reihe beriicksichtigt. Tut man dies bis zur zweiten Potenz, so erhalt man 


c _ V2gD [ku vo" — g D2 n) (1 —n)] + [g D(2 —n)P? —g D(2 —n) 
vo kw" vo” —g D (2 —n)(1 —n) F 
Am besten rechnet man den angenadherten Wert von c/v aus (13c) aus und setzt ihn auf den rech- 
ten Seiten von (13a) oder (13b) ein. 
Ist es nicht statthaft, Ff, zu vernachlassigen, so eliminieren wir dP/ds aus den Gleichungen (6) 
und (7), deren rechte Seiten zu null werden, und erhalten 
aoe v—c\?—n o4 DAA yee YU ee 
VslyeD° = = lee Spowil 1 er we : 


Daraus 1a6t sich entsprechend (13a) und (13b) eine Gleichung dritten Grades von der Form 


(13a) 


c 
eZ 
wv 
i 


(13c) 


GE: ie m= +p=0 (14) 
ableiten, worin 1, m und p anstelle der folgenden Ausdriicke gesetzt sind: im Falle 0,5 >n > 0 
pe 2g D (2 vg —# 71) Rat ies 8) ie A re 
Ye [hk wi” v” (1 —e/v)” —2g D] 4 Ly, [kwh v” — 2g D(1 —e/v)"] 
—2gDun 


Us yg [k wi” vo” (1 — c/v)” — 2g D) 


1 W. Barth, Fortschritte der Verfahrenstechnik (1952/53) S. 63. 
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und im Falle 1 > n> 0,5 


2gD e\i=n ,_ bg DLUn—y,_) 1 — ey) = Day wy 
SS 5 ar Gaon ( i aaa 4A Lv" wy” Vg k 


a 9 
kvm yet 


2%442gD (1 a 


n 2— 
ygkv wy v 


P = 


Die Koeffizienten mit Strich treten an die Stelle derer ohne Strich in (14). Wendet man wieder 
die Zerlegung in eine binomische Reihe an, die man hinter dem quadratischen Gliede von c/v 
abbricht, so ergibt sich ebenfalls eine Gleichung von der Form (14), deren Koeffizienten mit dem 
Index b versehen seien: 

gD(2—n) [27g —-e yn (l—n)] 


Ine Vg Lk wee v” —g D(2—3n-+ n?)] ” 
_ 82 DL [pian —7 ia Avitn we? v4 
Has AL y, [kwi-” vo — gD(2—3n + n®)| : 
= 2g Dun 
Py ~~), [k wi 0” — g D(2—3n +n] ° 


Oft wird man die Reibung zwischen Luft und Rohr gegen die zwischen festem Gut und Rohr 
vernachlassigen kénnen, also 2 = 0 setzen. Aus den durch Reihenentwicklung erhaltenen Ko- 
effizienten |,, m, und p, ergeben sich nur ungenaue Werte fiir c/v, so daB eine Nachrechnung 
mit den Koeffizienten 1, m, p, bzw. I’, m’ und p’ nétig ist. Zur Lésung von (14) wird man sich 
meist der goniometrischen Methode bedienen miissen, die fiir das oben ausgerechnete Beispiel 
fiir Staub sich in folgender Form ergibt!: 


(1/3)? —I m/6 + 0,5 p (15) 
[(l/3)? —m/3] (0/3)? — m/3 | * 


Fiihrt man hier die Koeffizienten J, m, und py fiir den Staub ein, so wird c/v = 0,864. Setzt 
man diesen Wert in die Koeffizienten l’, m’ und p’ ein, so erhalt man das genaue c/v zu fast 1. 


: 
ee cos | 60° + tare cos 


5. Die Aufgabestelle. Wenn man nach den oben aufgestellten integrierten Formeln von Anfang 
an die Zustande in einer Férderleitung berechnen will, fehlt noch der Ausgangswert cy. Vor allem 
deswegen seien die Verhaltnisse an der Aufgabestelle erértert. Sie sei so ausgefiihrt, wie das 
Schema in Abb. 3 zeigt. Die Anfangsgeschwindigkeit 
der Kérner in der Stromrichtung sei beim Einfallen co. 
Der Abstand eines Punktes vom Anfang der Aufgabe- 
stelle in Stromrichtung werde mit b [m] bezeichnet. 
Die gesamte Lange der Aufgabestelle sei b,. Die Zu- 
gabe des festen Gutes soll gleichmafig iiber b, verteilt 
sein. 

Die Kornchen haben also in einem Querschnitt 
verschiedene Geschwindigkeiten, so da®B Zusammen- 
st6Be vorkommen. Ich nehme an, daB alle zuerst dem 
Luftstrom zugefiihrten Kérnchen auf die spater dazu- 
kommenden stoBen und daB das feste Gut unelastisch 
ist. Unelastische Kérper gleichen dabei ihre Ge- Abb. 3. Schema der Aufgabestelle. 
schwindigkeit aneinander an, so daB ich voraussetzen 
kann, daB sich fiir jeden Querschnitt eine einheitliche Geschwindigkeit c’ ergibt. 

Ks erfolgen dabei auch schiefe StéBe, durch die Korner gegen die Rohrwand geworfen werden. 
Die Wand sei vollkommen elastisch. Dann tritt durch den Aufprall kein Energieverlust ein. 
Wiirde ein von der Wand zuriickgeworfenes Korn ohne andere Behinderung wieder mit dem 
zusammentreffen, von dem es den vorhergehenden schiefen Stof erhielt, und ware dieses von 
der anderen Seite der Rohrwand ebenfalls ohne Behinderung reflektiert, so wiirden beide K érner 
in Richtung der Rohrachse mit gleicher Geschwindigkeit weiterfliegen, wenn sie vor dem ersten 
Zusammenstof beide diese Richtung hatten. Die Gesamtwirkung aller schiefen StéBe wird der 
gleich sein, die erzielt wiirde, waren die StéBe alle gerade gewesen. Die Wand wirft die Kérnchen 
mit derselben Wucht wieder in den Strom, mit welcher sie auftrafen. In Wirklichkeit gibt es 


1 Siehe Hiitte 1948, I. Band, S. 77. 


XXIV. Band 1956 Hesse: Strémung in Blasleitungen 305 


weder ganz unelastische, noch vollkommen elastische K6rper, und man kann daher wohl an- 
nehmen, dafs sich der Fehler, den man durch die Annahme der elastischen Rohrwand macht, 
ausgleicht gegen den, der sich aus der Annahme der unelastischen Kérner ergibt. 

Es kommen aber auch exzentrische StéBe vor, durch die Kérner auBerdem in Drehung ver- 
setzt werden. Eine Nachrechnung ergibt, daf selbst bei Maiskérnern die dafiir verbrauchte Ener- 
gie sehr gering gegen die fiir die translatorische Bewegung ist. Die Drehung soll vernachlassigt 
werden. 

Die friher aufgestellten Gleichungen sollen nun fiir die Aufgabestelle abgedndert werden. 
Anstelle des gleichbleibenden G tritt das verinderliche Gutgewicht G b/b,. Die Stetigkeits- 
gleichung lautet also hier 

LE, Gb 


yi v | b. Hck seat F, Pe (I ) 


Wahrend von den Reibungsgleichungen (2) und (3) die fiir die Luft unveriindert gilt — anstelle 
von ds tritt db—, lautet die fiir das feste Gut 


ee Gb 
Pas koeDli Fe) abs (3b) 


Zur Aufstellung der neuen Impulsgleichung schneidet man eine unendlich diinne Querschnitts- 
scheibe von der Dicke db aus dem Rohr der Aufgabestelle. Zwischen Vorder- und Riickseite der 
Scheibe tritt das Gutgewicht G db/b, mit der Geschwindigkeit cy zu dem bereits vorhandenen 
G b/b, hinzu. Die Geschwindigkeit des neu hinzutretenden wird durch StoB nach Voraussetzung 
auf c’ erhéht, die des bereits vorhandenen durch die Férderluft um de’. Die Gleichung (4) wandelt 
sich also in 


£ Cie eC re tab Nyhan ah 
FdP — g dv 6, de z (c Co) b, | (Fa ve Yi+ kc c7] 2¢ D (4b) 
oder 
is G ; ; ‘ oni in 
F dP == dv + [d(b¢) — 0 db] +(FAe ! ke Gr) agp: 


Aus den beiden Impulsgleichungen fiir das feste Gut einerseits und die Luft andrerseits erhalt 
man durch Division mit dt = db/c’ die Kraftegleichungen 


Chdh Ch ween Cel tdbe— . b cf? 
Si pled EC, 61 
be yy db 5, ( wy || db ‘ aE Wnt 22D (6b) 
LRA Ee OU URC a Te dv, vy 
ee st F ie 71 
db bs a wy ot ap (a) 


Soweit in das Rohr festes Gut mit der Geschwindigkeit c) eingefiihrt wird, wird auBer dem 
Reibungsverlust zwischen den strémenden Stoffen, der sich entsprechend der Gleichung (8) 
hier zu ra 

R= (71) = 65 ST ob (8b) 

Re Dec wi 

ergibt, auch der Arbeitsverlust infolge der Zusammenstéfe der Kornchen S [mkg/s] in Warme 
verwandelt. Um ihn zu bestimmen, bildet man aus (6b) durch Multiplikation mit db die Arbeits- 
gleichung und zerlegt das Differential d (b c’) in zwei Summanden 


Gb Gb ft e\ 2 (xy, he G e fae db Gbhc? 
Eee ( = Deeg Oe ae og th, (16) 


Der erste Summand auf der rechten Seite der Gleichung stellt die Zunahme der Bewegungsenergie 
des vom Anfang bis b bereits aufgegebenen festen Gutes auf dem Wege db dar, der zweite Sum- 
mand die aufzuwendende Sto arbeit, die die Geschwindigkeit der auf dem Wege db neu hinzu- 


: : G¥o2 db 
tretenden Menge (G db/b,) von co auf c’ erhéht. Vor dem Sto war ihre Wucht S 7 ine danach 
ie ; ; G [c”? Co" 62% =— 2 cc, ¢n"\ db : 
a a Bringt man jenen Summanden in die form ; (3 5 as oat ae » SO er- 
kennt man, da der StoBverlust 
eu G (c’ —cy)* db (17) 


g Z b 


é 
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ist. Dies ist derselbe Betrag, den man erhielte, wenn man dS nach der Lehre der Physik fiir 
den unelastischen StoB berechnete. Ist c) = 0, so wird die Halfte des Arbeitsaufwandes fiir die 
St6Be in Warme verwandelt. 

Eliminiert man dP/db aus den Kraftegleichungen (6b) und (7b) und setzt man die dafiir 
erhaltenen Ausdriicke einander gleich, so erhalt man die allgemeine Beziehung zwischen c’, v 
und b 


ZC 


ee nme + 1)2g(S A) 42 bed tae es 


.. b dh bee we v, db ALY, D 
Ist das Volumen der festen Bestandteile klein gegen das Luftvolumen und nimmt man v konstant 
: GodP . Bs F 
an, 80 verschwindet ee (6b), und man erhalt unmittelbar 
e?g 
dc?) , ¢/2 c’ ¢ oy 64 \\2 TC 
db FBS 258 2 6( oa Das eo) 


Wenn n = 0 ist, laBt sich diese Gleichung zum Teil integrieren, indem man c’ durch db/dt ersetzt. 


Es wird unter Verwendung der Abkiirzungen a = Pe und h = g/w¢ — k/(2D) 


b 
Pee. 2a 1 ehb db 1 
ee L eal a) Pave [bets (19) 


0 


Das ungeléste Integral verschwindet fiir den Anfang der Aufgabestelle. Dort kann man 
de’ /db = c'/b setzen und bei ganz kleinen Werten von b aus (18) c’ berechnen. 
Wenn F, gegeniiber F klein und v konstant ist, gewinnt man aus der Impulsgleichung (4b) 


oe , 
aes WN oe 
0 


b 
Gb , 2 
AP aia eae oe AS few. (20) 
Darin ist nach (19) 


b b b 
. , 2 cf /ehd f 
[e d(b2) oa —5) ! 7 P(E ! i) a a] ! 2a | ar(femat ania 
0 0 0 


fir n = 0. Fir b = 0 wird das ungeldéste Integral wieder gleich null, und die Integrationskon- 
stante K berechnet sich aus diesem Fall zu 


Damit wird 


b "4 b 
; ee ohe = | ae ees 2 hb] 
le d(b?) ~ al — i) ! mut ! | ao ! L2av [on(frmoe)a. 
0 ) 0 


0 


Der Reibungsverlust innerhalb des Gemisches ergibt sich aus (8b) zu 
b b 


eh aa | ef (bt) = fe w) (21) 


0 


ab 
v3 lee. 


und der Sto®verlust aus (17) zu 


b b 
: G / ! , / 
Sm ze [o—26 | a+ ct a 
0 0 


1 Es ist d(c’*)/db = 2 d?b/dt®?. Die dafiir erhaltene Gleichung integrier 
der Variation der Konstanten und erhalt iach siihccumicde s  S SIE 


a 


ead 1 % av 
‘ “55 (8+ F| i om be ae a 


Mies 


Mee Ee die Integrationskonstante bedeutet. Nach der Anfangsbedingung ist ¢’ = cj, und { dt = 0 fir b = 0 
ultipliziert m i it i s “a 
pliziert man die Gleichung mit 6, so kann man daraus C = c//h — 2a/h2 bestimmen. 
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In Abb. 4 sind die Gutgeschwindigkeit c’ und der Druckverlauf P, berechnet nach (19) und (20), 

fiir Weizen bei denselben Annahmen wie fiir Abb. 1 fiir eine Lange der Aufgabestelle b, = 0,3 m 

dargestellt. ; , 
Die fiir das Ende der Aufgabestelle erhaltenen Werte sind die Ausgangswerte c, usw. fiir die 


anschlieBende Rohrstrecke. 
Die vorstehenden Ausfithrungen sollen da von Nutzen sein, wo man fiir eine Foérderleitung 
den Leistungsbedarf und die Abmessungen berechnen will bei gegebener Menge von Férdergut 


ie ie 


S «700d kg/m’ abs & 


iy 


12,9 


Weg b 
0 10 20 fem] 30 
bg 
Abb. 4. Zunahme der Gutgeschwindigkeit c’ und Druckverlauf P wiahrend der Aufgabe von Weizen ins Rohr bei denselben Annahmen 
wie fiir Abb. 1. Die Aufgabestelle sei ausgefiihrt, wie in Abb. 3 dargestellt, und habe eine Linge von b, = 0,3 m. Zum Vergleich ist 
die Gutgeschwindigkeit c aus Abb. 1 gestrichelt ecingezeichnet. oH = 0. 


und Luft. AuBerdem geben sie bei der Messung und Regelung der Gutbeladung durch den Druck- 
abfall einen Anhalt dafiir, wie lang die dazu nétige verengte Rohrstrecke sein muB. Wie ver- 
schieden diese Linge sein kann, ersieht man aus einem Vergleich der beiden Schaubilder fir 
Weizen und Kohlenstaub. Dieser fiihrt zu dem Gedanken, durch eine solche Messung des Druck- 
unterschiedes an einer verengten kurzen Rohrstrecke, die von einem Gemisch von Gas und 
Flissigkeitstrépfchen durchstrémt wird, die GréBe der Trépfchen ungefahr angeben zu kénnen, 


wenn die Mengen von Gas und Flissigkeit gegeben sind. 


(Eingegangen am 1. Februar 1956.) 
Anschrift des Verfassers: Dr.-Ing. Horst Hesse, Kéln-Flittard, Leverkusener StraBe 12211. 
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Einflu8flichen fiir die am bogenférmigen Rand eingespannte 
und am geraden Rand freie Halbkreisplatte 


Von H. Miiggenburg 


1. Einleitung. Die Halbkreisplatte ist eine von den Plattenformen, die einer Berechnung auch 
bei verschieden gelagerten Randern mit verhaltnismaBig geringem Aufwand zuganglich ist. Die 


Berechnung vereinfacht sich fiir jene Stiitzungsarten des geraden Randes, fiir die man Funk- | 


tionen finden kann, die die Bedingungen am geraden Rand schon erfiillen. Das ist u. a. dann 
miglich, wenn der gerade Rand frei ist, welcher Fall hier behandelt wird. : d 
Es sei vorausgesetzt, daB die Platte isotrop und homogen, die Plattenstarke h klein gegeniiber 
den iibrigen Abmessungen, die Durchbiegung w klein gegeniiber der Plattenstarke und das 
Hookesche Gesetz giiltig sei. Wenn E den Elastizitatsmodul und « die Querdehnungszahl be- 
deuten, ist 
Eh? 


a (i — 22) 


die Plattensteifigkeit. Die elastische Flache w wird durch die Differentialgleichung 


AAw(x, y) = ey 


bestimmt, in der p(x, y) die Flachenlast bedeutet.. Wirkt nur die Einzellast P = 1 an der Stelle 
Xo» Yor dem sogenannten Aufpunkt (Abb. 1), ist also p(x, y) = 0, so sei G(x, Yo % ¥) die zur 


ADP 


= =a 


Abb. 1. Abb. 2. 


Einzellast gehérige Biegeflache. Nach dem Maxwellschen Reziprozitatssatz ist die von der an 
der Stelle x, yg wirkenden Einzellast P = 1 hervorgerufene Biegefliche die BiegeeinfluBflache 
fiir den Punkt x, yy und zwar fiir die gleichen Randbedingungen. 


Die Funktion G(x, yo; x, vy) setzt sich zusammen aus einem im Aufpunkt singuldéren An-— 


teil G) und einem regularen Anteil G,, der die Randbedingungen der Platte zu erfiillen hat!: 
Oxo» ¥os XY) = Gol Xo» Yos % ¥) + Ci(%o» Yos XY) « (1) 


Die EinfiuB8flachen fiir die Schnittgré®en erhalt man, wenn man die BiegeeinfluBflache nach 
den Aufpunktskoordinaten xp, yy differentiiert. Der Regularanteil, der ein Integral der Differen- 
tialgleichung AAG, = 0 sein muf, hat die gleichen Randbedingungen zu erfiillen wie die Biege- 
einfluBflache?. 

Bei der vorliegenden Platte mit einem freien Rand werden die Regularteile so gewahlt, daB 
sie von vornherein die Bedingungen des freien Randes erfiillen, daB also das Moment 


F 02G 0G 
a Nl b 532 | = 9 


(2) 


1 Ph. Frank und R.v. Mises, Differentialgleichungen der Ph 


P. sik, Bd.I, S. 631ff., B hweig 1925. 
A, Pucher, Ing.-Arch. 12 (1941) S. 76. : ee 


_ XXIV. Band 1956 Miiggenburg: EinfluBflachen fiir die freie Halbkreisplatte 309 


und die Stitzkraft 


Ce (pet | = 0 
Z Ox® PY Ox by? | (8) 


am Rand cite 0 sind, wenn die y-Achse den freien Rand bildet (Abb. 2). Auch die Singularteile 
erfillen die Bedingungen des freien Randes, da sie mit (2) und (3) hergeleitet werden!. Die 
EinfluBflachen werden fiir die Querkontraktionszahl fe = 0,3 berechnet. 


4 2: Biegeeinflubfliche. Der Aufpunkt liegt in der Mitte des geraden freien Randes, der Koor- 

_ dinatenursprung im Aufpunkt (Abb. 2). Die BiegeeinfluBfliche ist symmetrisch zur x-Achse. 

Deshalb braucht nur eine Plattenhalfte betrachtet zu werden. Der Singulirteil lautet in Polar- 
koordinaten (Abb. 2)? 

“9 Ts eer’ USS oe ; r 1 

G(r. @) Seeman - | orl sin 2 g— cos 2. n> Pars; ; 


Fuhrt man mit dem Plattenhalbmesser a die dimensionslose Variable 9 = r/a ein, so wird 


(4) 


a2 ( 


a l+u ; 1 
Colo.) =7E7 yan |e Li“ (p sin 2 y — e082 yng) cea (5) 


Die Beziehungen zu den kartesischen Koordinaten lauten x = r cos Y, Y =rsing. 
Der Regularteil G, wird gebildet aus einer Summe von biharmonischen Funktionen v,: 


G, = PS Cn Un » (6) 


wo c, Freiwerte sind. Da die Funktionen v, schon den Bedingungen des freien Randes geniigen 
sollen, braucht nur noch der eingespannte Rand erfiillt zu werden. Das geschieht punktweise 
mittels der Konstanten c,. In diskreten Punkten miissen Wert und Neigung der EinfluBflache 
Null sein: 


ele en0) and Ese (7) 


Diese Bedingungen werden hier erfiillt fiir gp = 0°, 18°, 36°, 54°, 72° und 90°. Als weitere Be- 
dingung wird fiir y = 90° 

0G 
ee , 8) 


eingefiihrt, womit der Randverlauf der EinfluBflache in diesem Punkt eine horizontale Tangente 
erhalt. Dadurch werden die Bedingungen an den Ecken der Zweieckplatte streng erfiillt. Es 
wurde bereits darauf hingewiesen, daB dies wesentlich ist?. Dain jedem der vorgesehenen Punkte 
zwei Randbedingungen streng erfillt werden und fiir die Eckpunkte je die zusatzliche Be- 
dingung (8) bestehen soll, ist ein lineares Gleichungssystem mit 13 Unbekannten zu lésen. 

Die Funktionen v” werden so gebildet, daB sie zur x-Achse symmetrisch sind und den Be- 
dingungen 


0? v, 07 vp 

age ape |= 0 (9) 
Cou Our 
E t (2—yp) ax 3i3| fore 0 (10) 


geniigen. Sie werden zusammengesetzt aus je einem harmonischen Polynom P, und einem 
biharmonischem Polynom B,,. 
Man erhalt lee ae, (11) 


wo z, einen Freiwert bedeutet. Die harmonischen Funktionen P,, sind Realteile der komplexen 
Funktionen 

C (ae uy) 
wahrend die biharmonischen Funktionen B, die sogenannten normierten Polynome sind’. Der 
Faktor z, wird so bestimmt, daB v, fiir geradesn der Gleichung (9) und fiix ungerades n der 
Gleichung (10) geniigt. 


1 4, Scherer, Allgem. Herleitung singularer Losungen der biharmonischen Gleichung, sowie Aufstellung 
Greenscher Funktionen am Beispiel einer Dreiecksplatte mit freiem Rand, Diss. Aachen 1955. 

2 A, Néddai, Die elastischen Platten, Berlin 1925. 

3 F. Schultz-Grunow, Z. angew. Math. Mech. 33 (1953) S. 227. ; 

4 K. Zweiling, Grundlage einer Theorie der biharmonischen Polynome, Berlin 1952. 
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Fiir gerades n wird nach Einsetzen in (9) 


alee on [age tH cele | 

use | (12) 

Zn 2 2 
eer _ 

erhalten und fiir ungerades n nach Einsetzen in (10) 
a ee 

ee pte eae ae 

ly = 
Peie-o el, , 


Dabei erfillt dann v, bei geradem n auch die Gleichung (10) und entsprechend bei ungeradem n 
auch (9), wie folgendes Beispiel zeigt: v, wird zusammengesetzt aus P, = x* — 6 x? y? + y* und 
B, = x*— 3 x? y?. Unter Beriicksichtigung von (12) wird 


— 12y? (1 — 
ty = 20 -W. (14) 
Man erhalt so 
vg = (x4 —6 at 9? + y*) —2 (1p) (xt 3 ay?) = — (12) ot — 6 we ay? $y". (15) 
Setzt man nun y in (10) ein, so ergibt sich, wie behauptet, 
ha) Ov A 
aa + 2—#) a53,8 ee (16) 
Die Funktionen v, lassen sich komplex ausdriicken mit €” = (x +7y)". Fir gerades n gilt 
Un = Re 6" —= (1—y) Re€ Re tr (17) 
oder in Polarkoordinaten mit dem dimensionslosen Halbmesser 0 
tm = Q" [eos ng —-—> (lL—p) COS @ COs (»—1)9]. (18) 
Fir ungerades n hat man 
+4 n 
v, =>" Rel ++ (1—p) Ret Refr—! (19) | 
bzw 
n{t+u n 
et 10 5 cosny + (1x) cos g cos (n —1) 9). (20) 
Im Einzelnen wird 
ye 
v1 =O cosg, 
V, = @* [cos 2 p — (1 — p) cos? g], (21) 


il 
V3 =e" aH cos 3p ++ (1 —pi) cos p cos 2 9) ; 


v, = e* [cos 4p — 2 (1 — yz) cosy cos3 y] usw. 


Die EinfluBflache bekommt nun die Form 


0, 0; = e 2 ayy os 1 ’ 
G (0, 0; 0, ¢) TOF aN Ing + T— pl? Sin 2p — cos 2 p Ing) eer + Di enrn]. (22) 
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Da 13 Bedingungen zu befriedigen sind, lauft in (22) n von 0 bis n=12. Die Auflésung der 
13 Bedingungsgleichungen liefert 


C= 1,7876274, c, =  0,0129321, 
c, = — 3,9866612, cs =  0,0020971, 
c=  2,8091287, cy = — 0,0128670, 
cs—=  0,4260921, C9 = — 0,0078559, 
c,= 0,1942653, e,—=  0,0101417, 
c; = — 0,0503123, Co = —0,0103238. 
Cs = — 0,0377284, 


Nun ist die BiegeeinfluBfliche, die in der Abb.3 im Schragbild und in der Abb.4 durch einen 
_ Hohenschichtplan dargestellt ist, in dem die Héhen mit 2(3 +) N/a? vervielfacht sind, durch 
folgenden Ausdruck gegeben: 
eet WEN E (0,05 0,9) =¢'|Ing 4 1+# (9 sin 2 @—cos 29 Ing) + a 4 1,7876274 


l—p 
— 3,9866612 v, + 2,8091287 v, + 0,4260921 v, (23) 
+ 0,1942653 v, — 0,0503123 v, -—— 0,0377284 v¢ 


+ 0,0129321 v, + 0,0020971 v, — 0,0128670 v, 
— 0,0078559 v,, + 0,0101417 v,, + 0,0103238 r, . 


In Abb. 4 liest man den gréften Wert 1,7876274 a?/2(3 + )aN ab. 


ss 
a 


G05 Of 02 «=O% 06 08 10 12 Te 16 1787627 16 7% 2 10 08 O6 oF 


02 07 O05 
Abb. 4. 


3. EinfluBflache des Biegemomentes m,. Die EinfluBflache des Biegemomentes my, im Schnitt 
y = konst. wird ebenfalls fiir den Aufpunkt 0,0 ermittelt (Abb. 2). Sie hat die gleiche Symmetrie- 
und die gleichen Randbedingungen wie die Biege-EinfluBfliche. Als Regularteil G, kann daher 
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dieselbe Funktionenreihe 


GC, = DG, o, 
n 


(24) 


wie bei der BiegeeinfluBflache verwendet werden. Die Singularitat fiir das Biegemoment my 


fiir den Schnitt y = konst. am freien Rand lautet 
1l+yp TaN bed as 
Moy (3+ p)x Pelee = Camis 
oder in Polarkoordinaten (Abb. 2) 


eee rae 
Mo y = ere a — (1 —p) cos? o): 


a 


Man erhalt (25), wenn man den Singularteil der BiegeeinfluBflache 


= ae 1+ r 
Go (o> Yous % Y) a 216 Tpan™ _ (leo + %)?] In — 


ay 1 
—2 (x +35) (y—yo)are tg 1 at 7 


Xo 


mit r?2 =(« + x)? +(y—y)? fiir den Aufpunkt x9, yo in die Differentialbeziehung 


OG, 2G, 
oo ae n( aye ccs 
einsetzt. Da am freien Rand das Moment m), verschwinden mu (Abb. 2), gilt 
2G, 0G, 
me =—N (FE 7 ay 
oder = iS 
0G, eG, 
on Oye 
Setzt man (30) in (28) ein, so erhalt man? 
2G, 
UL = (1 ee) Bye g 


Die MomenteneinfluBflache m, formuliert sich somit zu 


1+ L r \2 n=12 
my (0, 0; 0,.~) = CEarAlte eI Sl) cost | a 2a Gn vf 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


Die Freiwerte a, werden so bestimmt, daB die Bedingungen des eingespannten Randes 


= omy me 
iene) aa ee] =o 


in denselben Punkten wie bei der BiegeeinfluBflache und die Bedingung 


Omy an 
a e=l a. 


(33) 


(34) 


fiir py =90° erfiillt werden, so daB wieder ein Gleichungssystem mit 13 Gleichungen zu lésen ist. 


Die Konstanten a, ergeben sich zu 


a = — 0,8688154, a, = — 0,0022810, 
d,= 0,7351904, ag = — 0,0030799, 
a, = — 2,0683212, a, = 0,0061746, 
ds = — 0,3138062, a, = 0,0097565, 
a, = — 0,4661287, a,, = — 0,0025074, 
a,= 0,0214218, a1. = —- 0,0044466. 
a,= 0,0160941, 


1 Siehe FuBnote 1 von S. 309. 
2 Siehe FuBnote 1 von S. 309. 
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Die MomenteneinfluBfliche my, die in der Abb.5 
linien dargestellt ist, wobei die Héhen mit (3+ 


Form an: 


(3 + w)% 
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als Schragbild und in der Abb.6 durch Héhen- 
)a/(1 +4) vervielfacht sind, nimmt folgende 


Ce eae 


Neg) 
m,(0, 0; 0, ) = 'm G (1 — 12) cos? p | — 0,8688154 + 0,7351004 », 


— 2,0683212 v, + 0,3138062 v3 — 0,4661287 OF 


+ 0,0214218 v, + 0,0160941 vg — 0,0022810 v, 


(35) 


— 0,0030799 v, + 0,0061746 v, -+ 0,0097565 v1, 
— 0,0025074 v,, — 0,0044466 v,. 


\ 


Abb. 5. 


a -_—-—— , 
_— =I . 
< a 


GD) 


ed 
+0005 +025 +05 


+7 +2 +343 +2 41 +45 +025 +005 


Abb. 6. 


4. EinfluBflache der Querkraft g,. Der Singularteil q), der Quorkes(temntlatiche qy fiir den 
Aufpunkt 0, y, am freien Rand x =0 wird ermittclt mittels der Beziehung 


Da die Stiitzkraft 


Cx = 


am freien Rand oder 


1 Siehe FuBnote 1 von S. 309. 


AC, , OG, 36) 
10 = nee ” Bx Byol * 
0G, 2°, | 6 (37) 
| Oxe et Ox OV? |x=0 
POD ig aC) (38) 
Get: oa CaP) Oxy Oy? 
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ist, wird, wenn man (38) in (36) einsetzt, 
z eG) 8 
Me nla—n) | ~ 3+y)2 


ses ON ge) geste: (39) 


r2 ! ri 


mit r? = (y — yo)? + 2. 
In Polarkcordinaten mit 9 wird (Abb. 2) 
2 sin y sin p cos? p é 
= 1— 1) |. 40 


Die EinfluBflache der Querkraft q, wird ermittelt fiir den Aufpunkt 0,0 (Abb.2). Sie ist anti- 
metrisch zur x-Achse. Der Regularteil mu8 daher im Gegensatz zu den Biege- und Momenten- 
einfluBflachen aus Funktionen bestehen, die antimetrisch zur x-Achse sind. Diese Funktionen f, 


werden wieder zusammengesetzt aus je einer Bipotentialfunktion B, und einer Potentialfunk- 


tion P,, wobei die Potentialfunktionen P, Imaginarteile der komplexen Funktionen 


re 1 Ad 1 
ee etat 


und die Bipotentialfunktionen B, die normierten Polynome sind. Man erhalt 
fro = Pn + %n Bn, (41) 


wo Z, einen Freiwert bedeutet. Die Funktionen f, werden so aufgestellt, daB sie die Differential- 
beziehungen 


Of, Of, 
Fe IE 338 oo 0 (42) 
und 
Ofna copa bs 
0x8 a (2 — F) x fa) a 0 (43) 


erfiillen. Das wird erreicht, indem man fiir gerades n f, in (43) und fiir ungerades n in (42) ein- 
setzt. Man erhalt folgende Bestimmungsgleichungen fiir z,: 
Fir gerades n 


oe OP, — | OB, OB, 
ax3 ao et Ox ml onc [ Giri?) Ox f.| = ue 
oP, OP, 
“ Ox? + (2 — p) Ox el, (44) 
La === — == 
eB, eB, 4 
Ee ge 2) Ox | 
fiir ungerades n 
[ @P,, | P, | ne 0°B,, i 02B,, 
on @y? leno! "| axe | M Oy? |x=0 ee 


a2P,, a2P,, 
| Ox? gia dy? i 
oe 2F dy? }Jx—0 


(45) 


Bei geradem n erfiillen die Funktionen f,, ebenso die Bedin i i 

é nm gung (42), wie auch bei ungeradem n 
die Funktionen fn die Gleichung (43) befriedigen, was sich in ahnlicher Weise wie bei fe Biege- 
einfluBflache nachpriifen la8t. Nachdem z, ermittelt ist, wird nun mit ¢" = (x-+1y)" fiir gerades n 


1 ae 
fr =" Im Ce Rein — (46) 


bzw. in Polarkoordinaten mit 0 


fn = 0" 


PSS ih 
aa sin n @ +- 


i —p~ ‘ 5 
3 cos sin (n—1)@/, (47) 
fiir ungerades n 


1 mn SS ¢ s . 
Sn == Im ¢ Stetore "Re €Im ta (48) 
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bzw. 

(pene I ee l—uyu é 

n= Q" |—- sin ——,— cos sin (n—1)q). (49) 
Beispielsweise lauten die ersten vier Funktionen 

fi=esing, 

te = e (sin ¢ cos Y); 

1 oe — 
he = 08 (5 sin Vie ne y cos psin 29), (50) 


1 Z a 
Sa es f sin A + ~ 7H cos p sin 39) : 


zane 
4 -3-Y-C-8WIl +086 +443 +2 Af ADS UE) 
Abb. 8. 


Fiir die EinfluBflache q, hat man nun den Ausdruck 


Wie i ; (51) 


n=1 


Ie : 4 
dy (0, 95 0,9) = one [sin p + (1 — 4) sin @ cos* al| a 


Die Konstanten b,, werden so bestimmt, da® die Randbedingungen in denselben Punkten erfillt 
werden wie bei der Biege- und MomenteneinfluBflache. Nun befriedigt aber bereits q,(0, 0; 0, ~) 
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die Randbedingungen im Punkt g=0 (Abb.2), so da® nur noch 11 Konstante 6, ermittelt 
werden miissen. Die Auflésung der 11 Gleichungen ergibt folgende Werte b,: 


Pps LE CRORES b, = — 0,0642880, 
b, = — 0,6733905, b, = — 0,0166327, 
b, =  5,6692867, b, =  0,0538504, 
b, = — 0,6405446, big = — 0,0157249, | 
b, = 1,9347084, b,, = — 0,1416655. 
b, =  0,0319381, 


Die EinfluBflache q,, die die Abb.7 als Schraghild und die Abb.8 im Hoéhenlinienplan zeigt, in | 


dem die Héhen mit } a(3-+y)a vervielfacht sind, formuliert sich endgiiltig zu 


EEE IEE G,(0;,.03°05 7) = > [sin p + (1 —p) sin ¢ cos? gp] + 2,4747744 f, — 0,6733905 fy | 


+ 5,6692867 f, — 0,6405446 f, + 1,9347084 f; (52) 
+ 0,0319381 f, — 0,0642880 f, — 0,0166327 f, | 


+ 0,0538504 f, — 0,0157249 f,, — 0,1416655 f,, . 


Die Rechnung hat gezeigt, daB zwischen den Punkten, in denen die Randbedingungen streng 
erfiillt sind, die Abweichungen der EinfluBflachen vom Sollbetrag gering sind. So betragt z.B. 
bei der BiegeeinfluBflache die maximale Abweichung 1,71°% bezogen auf den gréBten Wert der | 
EinfluBfliche. 


Die Anregung zu dieser Arbeit gab mit Herr Prof. F. Schultz-Grunow. Thm méchte ich an |) 
dieser Stelle meinen Dank ausdriicken. 


(Eingegangen am 3. Februar 1956.) 
Institut fiir Mechanik der Techn. Hochschule Aachen. 
Anschrift des Verfassers: Dipl.-Ing. Hans Miiggenburg, Aachen, Techn. Hochschule, Institut fiir Mechanik. 
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Eigenvibrations of Curved Elastic Bars According 
to the Method of Internal Constraints“* 


By E. Volterra 


1. Introduction. In a previous paper! the following equations of motion were obtained for a 
curved bar of length L, of curvature c and radii of inertia Tye ty: 


071, 0? 
igat ~ Ode + (A +26) eo — (24 36) cr 6 Gey, 


A 2G 
(+ )r 7 Ox Ox 
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Cri Ge — A+ 26)Ay—AS2 Aw ++ 8G) eM trem 4+ 6c eB 
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pe Cr Sea Cs ye. 
OP ly Ov». 2 Ou, 
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(A + 26) ri, 3 F — Cuz —G a2 —(A + 3G)e ry a" = 0 ry opp? 
i Cae Olin 072, 0?y, i OPA, 
Coa TOG TG cas =O aa — Oe aE 
together with the following boundary conditions at x =0, x =L: 
OA, OV, _ 
las teas — omy OA, = 0, 
pe tre Ge tem +A, Cie, 
OA f) 
ear Fe te Ad Oy 0, 3) 
Ob, 
oe)? rai 
OV, 1 2 OA, | J 
a+ 26) M2 + (A426) re F* A426) ery +Ady +4 ps] dre = 0 


* The material in the paper is based on an investigation, conducted at the Rensselaer Polytechnic In- 
stitute, Troy, New York, under the sponsorship of the Office of Ordnance Research, Contract No. DA-30-115- 
ORD-709. Project No. 454.13. 

1 E. Volterra, Ing.-Arch. 23 (1955) S. 402. 
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and 
OA, av.) « Ou 
et eae] =o. [oe + © Hs] Oy = 0 
Ou Ov, 2) OA, Se) 
peo Hy] Site = 9, FE + ae +6 r85E] Om, =O. | 


The above equations and boundary conditions were derived on the assumption that the follo- | 


wing equation of constraints is verified during motion: 


V(x, Ys %t) = V(x,0) + y A(x, t) + 2 a(x, 0) (5) | 


and by the application of Hamilton Principle !. 


In equation (5) V represents the elastic displacement vector at a generic point P of the bar, 


y, A and uw are three vectors, a priori unknown, functions of the coordinate x and of the time t. | 


The bar is referred to a system of curvilinear righthanded coordinates (0, x, y, 2) in which 


the origin 0 coincides with a generic point on the central (baricentric) axis of the bar, x with the | 


central axis, while y and z with the principal axes of inertia of the cross-section of the bar, 7 
and G are respectively Lame’s constant and the modulus of rigidity of the material of the bar. 


Equations (1) and the boundary conditions (3) define the five components /,, Ay, Uz, Vx and Vy 


of the vectors A, ys and » (eigenvibrations in the plane of initial curvature of the bar), while, 


equations (2) and the boundary conditions (4) define the remaining four components /,, Ux, My 
and y, (eigenvibrations which occur outside the plane of initial curvature of the bar). 


In the case of a straight bar (c = 0), the above equations reduce to 
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él, A omy 
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* It should be pointed out that equations (1) and (2) appl i i i 
; pply to the general case in which the curvature c is 
a Elven ee ae variable x, eee that dR(x)/dx [R(x) being the radius of curvature of the central 
axis ol the bar] 1s a finite quantity. In fact in deriving the equations of motion, th ht icati 
Hamilton Principle, in the given order of approximation thevtrnis Se ee ee 
de(x) _ 1 dR(x) 
dx “R(x)? dx 


can be neglected. 
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The object of the present paper is to show how solutions of equations (1) and (2) which satisfy 
the boundary conditions (3) and (4) can be obtained directly once that the solutions of equa- 
tions (6), (7), (8), (9), (10), (11), (12) and (13) (Fundamental Solutions) are known. 


Assume as fundamental solutions i. e. as solutions for the straight bar case the solutions 


i 
t ial!) ¢ 


ee iw(l) we 
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ia 2) 4 
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= das] 
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While assume as solutions of the equations of motion (1) and (2) the following expressions: 


Wee (Law + Sor La) Se 


iol!) t (18) 


and 


k=1 ry 
U ACN, 
Pa) = 3S (Mya + Sy Myla * 
k=1 kr / f (19) 
(I 
Uz(%; t) = BS (Ma Bi 2; Or Mz,) e* * 
k=1 r#k 


Where the coefficients «,, B,, 7 and 6, are supposed to be small quantities compared with unity. 
The summation indicated by the .») sign is over all values of r between 1 and oo except k. 
rtk 

The problem of solving the equations of motion for a curved bar derived by the ,,Method 
of Internal Constraints“ is therefore reduced to the problem of finding the expressions for the 
numerical coefficients a,, B,, 7, 6, and for the frequencies wf), wf! (r, k = 1, 2,3, ...) in 
terms of the eigenfunctions L,,, Lyx, Lek, Mar, Myr, Mik, Nuk, Nyt: Nek and of the eigenvalues 
wi), c?), wo), wo) (k = 1, 2,3,...) of the corresponding equations of motion for the case of a 
straight bar. 


2. Eigenvibrations in the Plane of Initial Curvature of the Bar. On the assumption that the 
quantity (pure number) cr, is very small compared with unity so that (r,c)? can be neglected, 
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equations (1) can be written under the form ! 


(426) 2h _ 67,6 Be (24 + 36) ert Ge — rhe AGE — Com, 
=o} SB 20n0 se 
Cees (21436) cP + 2Acm+A+3C)cd,=0-Gy— 20 hose, 
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G ry ae Ady —(A FG) hee iy ee 
(hate @) en) Oe er ce ee | 


The eigenfunctions of the fundamental solutions L,, and Nyb3 Lyk N,; and M,;, satisfy the 
following differential equations, boundary conditions and orthogonality conditions?: 
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ip . . 4 “| 
he eigenfunctions L,, and Nyx; Lyx, Nx, and M,, must be chosen in such a way that not only 


the b iti iti 
eet ey conditions (21) and (24) but also the boundary conditions (3) of the general problem 


De o instance in discussing the problem of flexural vibrations of a curved bar simply suppor- 
ed, the eigenfunctions L,, and Ny» must satisfy the condition of simple support: at x = 0, 


W e Suppose also that in the same o Ppp. qu v x? > ms 
rder of a roximation the antities cv c“v c“v, Can be ne lected 
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Nee (displacement equal to zero), 
Ai 
a=) (bending moment equal to zero). 


In order to satisfy at the same time equations (24) and equations (3) one has to assume as associate 
eigenfunctions Lyx, Nx and M,; the ones which satisfy the following conditions, at « = 0, 
ax = L: 
dN x5 
dx 


Lips = Mux = |(h +26) FH! +A Lig +4 Mya] = 0. 


In the case instead in which the bar is built-in at both ends, the eigenfunctions L,, and Ny, must 
satisfy the conditions of fixity at x = 0, x = L 
NE 0 (displacement equal to zero), 


£,5,= 0 (slope equal to zero). 


The associate eigenfunctions L,;,, Nj, and M,,, in this case are the ones which satisfy the following 
boundary conditions at x = 0,x = L: 


Lyx = Mn = Nyx = 0. 


By substituting equations (18) into equation (1’), dropping out second order terms, in view of 
equations (20) and (23) one obtains 


dLy, dM, 
—orz [ol] (L- sta a Ler\— (24 +3 G)er? tk ey) ae Cane 


= 075 [OfMP (Dar + So Le tee Ge Cosy [oP ? 


rk 


(26) 
—o [ol (Ny + Zar Ny.) Ovens Oe 2d OM pee(3 G)e ts 
rZ 
Se al its 2 NT arc Up [ow], 
rZé#k 

, dL. 

er (OPP (Lys +2 Brlyn) +A +26) pee day: 
=—o 7; [oa]? (Ey + 2 Br Ly\ +2 orc Nyn [op , 
r# 
—T; 0 [wp]? (Mir + 2B, M.-) +Ae Ny, = — ry 0 [oy]? (Mex + & Br Mz) tere) 
\ ioe r+k 


dNy k 


dx 


— 9 fo??? (New + » Br Ner) —=(1 +26) c 
r#k } 
=-=0[a))? (Nee + » B, Ner| + 2o0rcjoW]? Len. 
rk 


By multiplying the first of equations (26) by L,, dx, the second by Ny; dx integrating between 0 
and L and adding them together, in view of the orthog onality conditions (22) one gets the value of 
@mWb 2 T? 
(OPN = wee ees 
given by equation (28). : 

Multiplying the first of equations (26) by Lz; dx, the second by N,, dx integrating between 0 
and L and adding them together one obtains, in view of the orthogonality conditions (22) the 
value for «, given by equation (29). ' 

Multiplying the first of equations (27) by Lyx, the second hy M,,,, the third by Nek integrating 
between 0 und L and adding them together, one obtains, in view of the orthogonality conditions 
(25) the value of . : 

Oop Ee 
[OO], = [M% a Q 
given by equation (30). 

Multiplying the first of equations (27) by Ly,, the second by Mas the third by N., inte- 
grating between 0 and L and adding them together, one obtains in view of the orthogonality 
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conditions (25) the value of B, given by equation (31): 
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In equation (28), (29), (30), and (31) 6 =r,/L. 

3. Eigenvibrations Outside the Plane of Initial Curvature of the Bar. 
the fundamental solutions L,;, and Myx; M,);, and N,,, satis 
boundary conditions and orthogonality equations: 

da? 
Gri — 6 Lin — G6 My, =— 0 F* [o®]? Lyn, 
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fy the following differential equations, 


(32) 
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As in the previous case, the ei i 

Sh > genfunctions L,,, My,, and M,;, N,, must satisfy not only the 
Pe conditions (33) and (36) of the fundamental cases but also the boundary conditions (4) 
yo = Saree problem. By substituting equations (19) into equations (2), by dropping out 
Second order terms, in view of equations (32) and (35), one obtains the following equations: 
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By multiplying the first of equations (38) by L.. dx, the second by M,, dx, integrating between 

0 and L, adding them together, one obtains in view of the orthogonality conditions (34), the 

following value for 


(TI) ]2 2 
(II)]2 —[%% Its) L? @ 
[Dy des Sa) E 


1 ‘ 
[ow], = z nee 


18+ | eNes Lena 
éj 


(40) 


ib 
2 d? Nuk 3—40 2 dM, 
é b 
By multiplying the first of equations (38) by L,, dx, the second by My, dx, integrating between 
0 and L, adding them together, one obtains in view of the orthogonality conditions (34), the follo- 
23 
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wing expression for y;: 
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By multiplying the first of equations (39) by M,;, dx, the second by N., dx, integrating between 
0 und L, adding them together, one obtains, in view of the orthogonality conditions (37), the 
following expression for 
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By multiplying the first of equations (39) by M,, dx, the second by N,, dx, integrating between 0 
and L, adding them together, one obtains, in view of the orthogonality conditions (37), the follo- 
wing expression for 6,: 


L 
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In equations (40), (41), (42) and (43) 6, =r,/L and 6, =r,/L. 


4, The Same Problem in Two-Dimensional Elasticity (Plane Strain). In two-dimensional 
elasticity (plane strain) on the assumption that the quantity (pure number) cr is very small 
compared with unity, so that (cr)? can be neglected, the equations of motion derived in a pre- 
vious paper? for a curved bar of length L, of constant radius of inertia r, and of curvature c can be 
written under the following form: 
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1 F. Volterra, 1. c. 
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with the boundary conditions at x = 0, x — L: 
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In the above equations /,, Ay; and yy are the components of the vectors 2 and? on the x and y 
axes. x coincides with the central (baricentric) curvilinear axis of the bar, while y with the normal 
to the central axis of the bar and is directed positively towards the center of curvature. E and o 
are respectively Young’s modulus and Poisson’s ratio for the material of the bar. 


Assume as solutions of equations (44) and (45) the following expressions: 
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r#k 


rales t) = B (Naw + 3 Br Nareien', slat) = 3 (Nya + 3a Nye} of, 
1 r¢#k 


k=1 rk 


(46) 


where L,,, Lyx, Nx~ and Ny» are the eigenfunctions of the fundamental case i. e. are defined by the 
same equations (44) and (45) when c = 0. 


The values a,, 8, are quantities which are supposed to be very small compared with unity 
and of the same order of magnitude of cr. 


The functions L,,, Ly,, Nx, and Ny; satisfy the following differential equations, boundary 
conditions and orthogonality equations: 
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By substituting the expressions (46) into the differential equations (44), by dropping out second 
order terms, in view of equations (47) and (50), one obtains the following equations: 


E(3 +0) dLy t, 
2 (1)]2 ce Oy Ler\— re | 
—or [a,”] Rte jee? Ales (53) 
E -s Bnei Ts I Foes + of2 P00 New| 
_-- 2(1 +0) Cc Nek = Or or ( xg Ib + 8 x ) H k C 


apy 


326 Volterra: Eigenvibrations of Curved Elastic Bars According etc. Ingenieur-Archiy 


E E(3+0) dNyy 
—o [al (Ns + 2 Or Ny) Se GS Ly, 4 2 (1 0) (ars 


Gone 


=—0 or (Nye + >) a, a) +2roca? I beep: 
r£k 


and 
: ; = E dL, Eo 
—or [o\?)]? (Eye + a pr Lye) ~ (eso) ree ae a = c Nyt 


=—or? ook (Lyk + aie bys) +2rocaz ieee 
rtk 


E dN. 
Q [wl?]? (ae Sp = fr Nar] <= F — ps = 0 on Net = yy Br Nar) 
T#k = 


r#k 


E 


SEA auc iyal Dyk 
By multiplying the first of equations (53) by L,; dx, the second by Ny, dx, integrating between 0 
and L, adding them together one obtains, in view of the orthogonality conditions (49), the follo- 
wing expression for 
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By multiplying the first of equations (53) by L,, dx, the second by N,, dx, integrating between 0 
and L, adding them together one obtains, in view of the orthogonality conditions (49), the follo- 
wing expression for a,: 
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By multiplying the first of equations (54) by L,, dx, the second by N,;, dx, integrating between 0 
and L, adding them together one obtains, in view of the orthogonality conditions (52), the follo- 
wing expression for 
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By multiplying the first of equations (54) by L,, dx, the second by Nx. dx, integrating between 0 
and L, adding them together, one obtains. in view of the orthogonality conditions (52), the follo- 
wing expression for f,: 
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5. Conclusions. In the present paper it has been shown how it is possible to discuss the 
problem of vibrations of curved bars once that the solution of the corresponding problem for the 
straight bar is known. The method of approach in deriving the equations of motion being com- 
pletely contained in the constraint equation (5) and in the application of Hamilton Principle, 
while the method of solving the equations being based on a perturbation method which utilizes 
as first approximation the solution of the straight case. 


More general dynamic problems of elongated solids in which the inertia characteristics of the 
cross-section of the bar are known functions of the variable x can be attacked and solved by the 
same method which has so far been demonstrated for the simple cases of bars having constant 
inertia characteristics. 


6. Appendix. In order to prove the orthogonality properties expressed by equations (22) 
or (37) in the eigenfunctions L,,, Ny, or Mzx, Nzx let 


Lyp = AE Mi ohh 
or 


De as py = 1 Ae Nay eee 


Then equations (20) or (35) can be written 


d?X,, 


a, 
dx? 


dx 


PY, | aXe 
dx? dx 


A+26)r COG =—ooir X,, G = —— OME Yme (09) 
Q 


Multiply the first of equations (59) by Xj dx, the second by Y;dx and integrate between 0 
and L. Adding them together, one obtains, in view of the boundary conditions (21) or (36) 


y L ‘i 
dX), dXj 
—wot([ 1 Yate + [x25 a) G4 20) yar raed | 
0 0 0 | (60) 
r j “ay, aY. ray 
"dY;, , | 4Y¥y 4Yj 2 i 
—6 [Xi Xjde—G | EEX; dx ON as cd di c | ig Xe dx | 
0 0 () 0 
In the same way one obtains for the jth mode 
L L is oe 
—eo3( [¥ Y,dx%—- re [XX l= —0 +26)r [eae sp dx | 
0 0 0 (61) 


0 0 


L L L L 
dYj dyad Ya dY, 3 | 
—6 [X)Xide—6 [ GP Xvdx—6 | Garde — 6 “ie Kide-| 
0 0 
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By subtracting equation (61) from equation (60), one gets: 
L L 
0 (wi — ok) ( Y, Yjdx +r? f DG 9.& as) = 0 
0 0 


which if @, 4 @; proves the orthogonality property expressed by equations (22) and (37). 

In order to prove the orthogonality property expressed by equation (25) in the eigenfunctions 
Lyk, Mz, Nuk, multiply the first of equations (23) by L,; dx, the second by M,; dx, the third 
by N,j; dx and integrate between 0 and L Adding them together, one obtains, in view of the 
boundary conditions (24) 


L 


2 2 
—— Ci QO [: | 
0 


L L 
Lyx Ly; dx +93 [ Me iV det | Nes ING; a 
0 0 


L L 
dLy7, aby; 
eee 2 | ee dae (1a 26) | Lyn Lyj dv 
d 0 
L L 
“dN, fe 
A [Spt hyjde—2 | Myp Ly; dx 
0 0 
; i - (62) 
— 615 [ETE A | yn Maj dx — (A +26) [ Mew May ax 
x x 
0 0 0 
r r aN dN s 
dN xt; ; xk xJ dN, j 
—A { ae M,jdx— (+26) | hale dx a Tn Ly, dx 
0 0 0 
L 
—} [eM dx 
dx 
‘ J} 
In the same way one obtains for the jth mode 
L L L 
— oho [ tet dx +7} | Ma M,1, dx +| NG a) 
0 0 0 
L L 
Mb poe GUL, 
Sse “de de A426) | Ly Lynde 
0 
* aN r M,; dM. 
xJ f d zJ Zz 
—1 | 2 es | M,; Ly, dx — 67 | re (63) 
i) 0 0 
L L TT; 
A | Ey Mex dx—(h +26) [ Maj Mex dx —a [ “7 Mu, dx 
Q 1) 0 
"aN; d r 
=— * dNxj INek dN xk Ne 
0 0 


By subtracting equation (63) from equation (62), one gets 
2 2 2 ” 2 7 ty 
QO (a; —— i) rf Lt Dy, dx ++ Tali M.; M.x dx +- if Ney Nuk as) = 0 
é 0 
which if Wj F Wp proves the orthogonality property expressed by equation (25). 


; In order to prove the orthogonality property expressed by equation (34) in the eigenfunc- 
tions Ly, and My, multiply the first of equations (32) by L,; dx, the second by My, ; dx and inte- 
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grate between 0 and L. Adding them together, one obtains, in view of the boundary condi- 
tions (33) 


L 


L L 
Laie. Lae Lej +1} Myx My;) dx = — Gri { Set Vai gy G [ bale, dx 
} | 
L L L L 
f dM, dMy; aii 
Pike M,,dx —G | M,, M,; dx —Gr ‘| ede c| Myx My; dx \ (64) 
0 0 0 
L 
2s [ Bae Myy ae. 
0 
In the same way one obtains for the jth mode 
L L 
: Ae db 
— of | (EL; L,, +71} My; My.) dx = —Gr? [a ri “* dx eis 
é ; 
r dM, ; dM 
mye [ te Mypde—6 | My; Max dx —6 15 i eae 6 | Mat dx | (65) 
0 6 ty 
L 


—6 | Lj My ae. 
0 


By subtracting equation (65) from equation (64), one gets: 
L 
@ (cof — oot) f (r2 Lex Ley +7} Myx My,) dx =0 
which, if w;  w,, proves the orthogonality property expressed by equation (34). 


(Eingegangen am 6. Februar 1956.) 
Anschrift des Verfassers: Prof. E. Volterra Troy, N. Y., 1567 Tibbits Ave., USA. 
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On the form of a jet issuing from a swirl atomizer 


By L. Valdenazzi 


1. Introduction. A swirl atomizer is a simple device for breaking a liquid into small droplets 
and consists, in its simpler form, in a cylindrical (or ,,vortex“‘) chamber with tangential inlets 
and a central outlet (fig. 1). The liquid is introduced into the chamber through the tangential 
inlets and issues from the central outlet with a tangential velocity component. This component, 
when sufficiently strong, causes the liquid to issue as a conical surface, formed by an unbroken 
liquid sheet (or ,,tulip“*) near the exit hole, and by droplets after a certain distance. The droplets 
are the result of the breaking of the liquid sheet under the 
action of the internal turbulence. 

The experience has shown that the axial length of the 
unbroken tulip is related to the velocities within the liquid 
motion, i. e. to the pressure drop through the atomizer. 
The higher is this pressure drop, the nearer to the outlet 
hole does the jet break into droplets. 

The liquid surface formed by droplets, in absence of 
external disturbing forces, such as gravity, the action of 
movements of the gas into which the jet is discharged, etc., 
is very nearly a circular cone, while in the case of unbroken 
tulip the surface tends to converge towards the jet axis; in 
the case of very small pressure drops through the atomizer, 
the tulip can actually become a bubble (fig. 2) without being 
able to break into drops. 


tangential inlets 


increasing pressure arop —=— 


Fig. 1. A swirl atomizer. Fig. 2. Form of the jet issuing from a swirl atomizer. 


A theory describing the form of the liquid cone has been put forth by Soehngen and Grigull }. 
In this theory the cone angle (¢ in fig. 1) is related to the ratio of the tangential and the axial 
velocity components at the atomizer outlet; both viscosity and surface tension of the liquid 
are not taken into account. The theory fails to agree with the observed results, even when the 
velocities are sufficiently high to support the assumption regarding the viscosity and surface ten- 
sion. The main reason is to be traced down to the fact that the Authors have evaluated the velo- 
city components at the outlet without taking into account the flow losses within the atomizer. 
When using correct values for the velocities, the Soehngen and Grigull theory gives better results, 
but only in the case of high pressure drops through the atomizer, i. e., when the jet breaks into 
droplets immediately at its exit. In the case of unbroken tulips, the influence of the surface ten- 
sion is no longer to be neglected. 

The present theory takes into account the surface tension, in addition to the liquid exit 
velocities, but not the viscosity. This latter, as the experience points out, does not need to be 
accounted for, because, if it were too high, it would not allow the very existence of a liquid 


' G. Scéhngen, Forschung auf dem Gebiet des Ing.-Wes., 17, 3 (Oa) pen ie 
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cone. In other words, the tulip forms at Reynolds number of the liquid flow sufficiently high 
to be neglected. 


2. Symbols. Coordinate system: cylindrical, the jet axis being the system axis, 


fi == radius, 
% == distance along the axis from the atomizer outlet, 
@ = angular coordinate, 
oe = 0147), 

s = liquid surface thickness, 
© = velocity (vector), v = velocity (scalar), 
oes | 

Vim 
E = total energy of the unit mass, 
M = mass, 
uN — s€e fiz. 5, 
Q =] I. = > 

he 
i Fre 
S = surface, 
V = volumetric rate of flow, 
Taga 

iy 

2 

We = Weber number — So 
eee Bee 11g. 2, 
€ = conical surface angle (see fig. 2), 

= 
Gaia oe 
@ = density, 
o’ = mass of the unit surface, 
o = surface tension (twice the experimental value), 
We see 118. 
Consistent units throughout. 


Subscripts: u = periferal, r = radial, m = meridional, z = axial, 1 = atomizer outlet 


(z= 0). 


3. Theory. a) Assumptions. The fundamental assumptions of the present theory are the 
following: 

a) The liquid motion is permanent; 

b) The liquid is inviscid; 

c) The fluid motion is a constant energy motion; 

d) The fluid into which the jet is discharged does not interfere with the motion, and it is 
at uniform pressure; 

e) The liquid motion is symmetrical about the axis; 

f) There is no velocity gradient along a radius of the revolution liquid surface, i. e. the velo- 
city is a function only of the distance from the atomizer outlet. 


b) The forces acting on the liquid. 1. Inertia forces: We will use a cylindrical coordinate 
system. Considering (fig. 3) an element of the unbroken liquid surface ABCD, let us investigate 


its equilibrium under the applied forces. 
The inertia forces have the following components: radial component 


dv, v2 
dM (a — oh (1) 


tangential component 
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axial component 
dv 
eZ 3 
dM ( i) 5 (3) 


dM being the mass of the surface element. The expressions in parentheses are the components, 
jn the chosen coordinate system, of the acceleration vector dv/dt. 


2. Surface tension forces: As external forces we must consider only those caused by the surface 
tension, as the pressure, according to the assumptions made, is constant throughout. The 
surface tension forces acting on the boundaries of the surface 


element are ! 
along AB 

PURI) Stet: le +(E) > 
along DC 

ein eee i aa 

—oa V (dz)? + (dr)? = — odz \/1 + i i 
along AD , 
aordg, 
along BC 
Fig. 3. An element of the liquid cone surface. =e; (r a dr) dp aT: do. 


Of such forces we need the components along the coordinate directions. 


As regards the forces acting on AB and DC, we may observe that they are tangential to 
the fluid surface as well as normal to AB and DC; they must therefore be tangential to the 
,-parallels“ of the revolution surface. Fig. 4 represents one of such parallels, with both resultants 
of the forces acting on AB and DC; we may deduce therefore that the component along the 
radius of these forces is 


2 
== 2 & Oks yi aL a cos (90° — a 3 


dr \? 
Seals Se in eS 
a \a,) ar. (4) 
2 
ar) C [ 
At LZ : SS 
\ / Z | 
/ 
ee 
ee 
\/ 
V 
Fig. 4. The forees acting along the parallels of the liquid Vig. 5. The fore:s acting along the meridians of the liquid 
cone surface. cone surface. 


The forces acting on BC and AD are tangential to the revolution surface as well as normal 
to its ,,parallels*‘; they must be therefore tangential to the ,,meridians“. The resultant of said 
forces is directed towards the curvature center 0 of the meridians (fig. 5) and has the value 


N =ordgy dy 


deducted as before. N may be decomposed into two forces, one of which acting along the radius 


Ncosa = WN 


5 OS Oe 
oe ie 


fs * The surface tension o must be taken twice the experimental value, the surfaces limiting the liquid sheet 
veIng two. 
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and the other acting along the direction of the axis 
dr 
; dz 
— Nsinao = —or dp dy — 


coma far\? 
a 
/ om dz 
We must now eliminate dy from the preceding expressions; as the curvature of the meridian, 
in the point considered is given by 


dr 
dz? 
|. /dr\2)3/2 
b+) 
as well as by 
as) ae 
dr\2 i 
dz y ++ (z:) 
we can deduce that 
d?r 
2 
dy = dz We 


d2r 

dz 

(i | 
and that along the direction of the axis 

dr d?r 

dz dz 


1 (eRe (6) 
b+ (a) 
Concluding, the components along the coordinate directions of all surface tension forces are 
along r 


d?r 


dr\? | 1 dz? 
—odzrdp yi + (&) - dr? 2? 
dz 


along z 
: dr d?r 


| (dry? dz dz? 
—odzrdp 1+ (i) Ate 
p+ (a) 


c) The equilibrium equations. If 0’ is the mass of the surface unit area, the mass dM of 
the element considered is 


d 2 
aM = 0' der dp|/1 +(3) 


the equilibrium equations may therefore be written in the following from: equilibrium along the 
radius 


| ar 
dv, v3 aj cel (7) 
Ht Sa eye n+(Z)i 

dz 
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equilibrium along the periphery 


dv, v, Vy 
dt ae 0, (8) 
equilibrium along the direction of the axis 
dr dr 
dv, Cam zm ze (9) 


Tin yes dr\?)2 
b +(z) | 
We must now observe that 9’ is not a constant quantity. In fact, as the surface thickness is 
given by 
y 


SS ioe 
V being the volumetric rate of flow out the atomizer, and v,, the meridional velocity: 
Om = Yor 8, 
the mass of the surface element considered is given by 
dM =osdS = oer dp ds yi el 
whence 


F oV M’ 


@ BLP Oa ho. 


As the motion is permanent, V is a constant quantity, but v,, il a function of z, so that also 
o’ is a function of z. 

The boundary conditions of the differential equilibrium equations are represented by the 
known values of the velocity components at the atomizer outlet, viz. 


Ur S= 18) 9 Pn) Vy = Vyu> Vz = Um = VUim> 


as well as the thickness of the liquid surface at the outlet. 

d) The integration of the equilibrium equations. 1. The elimination of the peripheral 
velocity component: Let us first consider eq. (8). As 

dr 


v, 


dt 


we immediately have 
dv v 


u u 
Ur dr + eee 0 


whence 
Vita g Tee (10) 
This result enable us to eliminate v, from eq. (7) and eq. (9). Introducing eq. (10) and M’ 
into eq. (7) and rearranging, we have 
d*r 
“2 o 1 dz? 


ve. Uses ae, 1 
Oe, Gi ee i M’ : r TPA uD 
; +( )| 
dz 


similarly eq. (9) becomes 


dr d?r 
v, du, o dz dz 
Vm dz a M’ i dr 212 : (12) 
b+ (z) 
dz 


Eq. (11) and eq. (12), when solved, would give v, and v, as functions of r and z. As the ulti- 
mate aim of this theory is the equation of the fluid surface, we need a third equation to eliminate 
v, and v,. This additional relationship is represented by the differential equation of the fluid 
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trajectories, viz. 

V, dr 

Ho de (13) 


2. The energy equation: The direct solution of our problem, i. e. the simultaneous solution 
of the three differential equations (11), (12) and (13), presents great analytical difficulties. 
Fortunately, our initial assumption of constant energy motion enables us to overcome such 


difficulties. 


Let us consider the energy stored in the surface element of fig. 3. There is a kinetic energy 


2 
os ds Ee 
and a ,,superficial‘ energy 
odS 
the sum of both being a constant, i. e. 
v2 ve 
z esdS > +odS =os, dS, > +od5,. (14) 


If we follow a given element during its motion, its volume remaining unaltered, we may write 
Sito = Sy 


so that, dividing eq. (14) by this quantity, and rearranging, we obtain 


v o> vi o 


wa Co | 


1 


The continuity equation ruling that 


we may finally write the equation energy in the form 


v2 o v? o 
yt apt em = a + ap tm 0?) 
that we shall use from now on. 


3. The solution of the equilibrium equations: Let us consider again eq: (11), (12) and (13), and 
eliminate the time t from them. As 


ss dr £. dz dm 
aide Geen apt fa 
and 
dm =a + (dr)? = ai +{ A 
we may write 
dr 
i dr “dz 


dz = 
Um > dm ; dr\2 
3 +(z) 
and substituting the latter equation into eq. (11) and eq. (12), we obtain the following equations in 
which the time variable is no longer present: 


Afr d?r 
dv, dz 2 aaa! dz* (16) 


pee Tr r) 
ae fat mt. Mere a 
é y+(z) E+ (& 
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du, 1 o dz dz : (17) 


We must observe now, as 
2 
2 2 2 % 
v = Um + Vu = Um + 


that the energy equation (15) represents a relationship between v,, and r. It may be used there- 
fore to eliminate one of the eq. (16) and (17). 
By differentiating the energy equation (15) we first obtain 


vdv + Ww (r dum + Um dr) = Vy Avy + Vm d¥m + ay (r du, + Um dr) = 0 


whence after some rearranging, we get 


x? o ‘ 
Nee 
doe eee (18) 
oO 
Ee Ge | 


Moreover, as the velocities are tangential to their trajectories, we may write 


1 
UV, =U, COS O =U, Fae 
i 
els 
y1+(@) 
dr 
, dz 
Vp = Um SNAG = SOS rae. 
Te 
1 hah 
y2 + (az) 
As 
dv,, __ dv,, dr 
dz dr dz’ 
the derivative of the first of these equations may be written in the following form: 
dr dr 
du, dz dv,, : dz 
dz dr\2 \ dr ee dr\? 
pre 1 — 
a (7) | i (z) | 
whence 
dr d2r 
dv, 1 dz dv, d2 


ay me PRC) dr\? | dr ™ dr\?}{ ° ™ 
pperey Bee Sees 
ye (i) P+ (a) 
Substituting now eq. (19) together with eq. (18) into eq. (17) we finally obtain, after some 
rearrangement 


d?r x2 o 
dz? mm 
ARNO o = f(r) (20) 
ie (za) rt ae 
dz M”2 


which is the differential equation of the liquid surface. Provided the function f(r) is known, 
its integration is quite easy. In fact, as it is 

d fdr\? _ 9 dr d?r 

dz \dz} ~~ dz dz?’ 
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eq. (20) becomes 


Integrating between r, and r we obtain 


in ae (z) | 49 | On (21) 


my 


whence 


i / ay f(r) dr 
ec 1 =) (22) 
and, after another integration 


r 


, a dr § 
sac ce 


Ty 


which is the required analytical equation of the liquid surface. 


4. The function f(r). We must now eliminate v,, from eq. (20) in order to obtain the function 
J(r). The relationship between v,, and r is given by the energy equation (15). If E is the constant 
value of the total energy 
v2 o 


| | 4 
a tg tb agtin = E 


we have from this equation 


o o Ree 
vy = MW’ | rat Vee r? & 2) . 


. : . . . > 2 . . 
The quantity between parentheses is certainly negative, as a < E. The radical is therefore 
if’ 


oO § : : ne hota : 
greater than w” and we must choose the + sign, if v,, is to be positive. Substituting v,, into 


eq. (20) we have 


S38 (24) 


Ee ‘fe? 2E os o = Bee =i 

vig (5 =) 2a Vga —(]-28) 

It is obvious that the expression of f(r) is too complicated to permit the effectuation of the 
integration of eq. (21) analytically. The solution of our problem must be necessarily numerical. 


5. Reduction to dimensionless quantities: By reducing the worked out relationships to dimen- 
sionless form, we have the advantage of bringing the number of the indipendent quantities to 
a minimum. 

A simple control of the dimensions of the various formulae shows that the integral 


[f@) dr 


is already a dimensionless quantity, as well as y(r). f(r) has the dimensions of the inverse of a 
length, and the second fraction of eq. (24) is therefore a dimensionless quantity. We want now 
to reduce to dimensionless form also the quantities entering into its expression. Observing that 


ee 2 Yim 8320 
2G 206 


MW’ 


= 0 Vim Sy Mh 


and assuming 
s&s, = UT 
i, a 
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where jz must be considered a known quantity, we obtain 


oO o 1 


a ot 
M Oita ety 


whence 
Oo 1 Vim 


We = Wee 
D 2 4 
the dimensionless quantity We being the ,,; Weber number 


2 
2 Vim ™ 
Lo = es 


o 
Let us write moreover 
aeRO fea R 
Dre r 
the expression of the energy E then becomes 
2 i ve 
vy o Vim vy u | lm 
Et yi al ! 5 Wee 
whence 
DAG: eek 
eee Se 
fir ae valle 


By introducing all these quantities into eq. (24), and rearranging, we obtain 


C? R 1 R? uy 2 | 
oe fet waa awa t Cl Rt) We | 
LOS 1 2 | 2 / R? moa 2 25 
Rlawat ( pC+l+—a\—ae wee ga) ti tay | Ga 
= —f'(R). | 
Integral (21) becomes 
r R 
dR ,, 
2[ fo) ar=2 [ Spe) (26) 
ra 1 
and writing 
zB 
Zi = Th 
we have 
ER 
dr dR / ad) 
= a mail = 
de dl : 
and finally 
R 
Z dR 
a RR 
\/ al Re F(R) 
1 el =, l (27) 


which is the dimensionless equation of the liquid surface. The parameters from which the func- 
tion Z = Z(R) depends are , C and the Weber number. 


6. The case of very great Weber numbers. When the Weber number is very great, i. e. when 
the surface tension forces are negligible in comparison with the inertia forces, we may write 


i) ee 1 
Fale 2) 
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and eq. (26) can be integrated analytically: 


ee ie gene 


— Q? 
where 
1 tl 
e == 1 *-- CG 
Substituting this result into eq. (27) we finally obtain 
/R?—1 
Vip a ae (28) 


In this case the meridional section of the liquid surface is a hyperbola (fig. 6). The inclination 
of its asymptote with respect to the axis is e. 


Eq. (28) has been found, as stated in the introduction, by Séhngen and Grigull in their 
elementary theory of the liquid surface form. 


/ 


R 
1 
ie ; 7 
Bs —— We increasing 
Fig.6. Form of the jet when the surface tension Fig. 7. Form of the jet at various Weber 
is negligible. numbres. 


4. Applications. As we have stated before, the determination of f’(R) is numerical. When 
the values of f’(R) for a given range of R values, starting from unity, are calculated, the inte- 
gration of eq. (26) is carried out numerically or graphically, and presents no difficulties, nor 
the subsequent calculations resulting in the Z values. 


For given values of y and C, it is interesting to analyze the influence of the Weber number 
on the jet form. A number of liquid surfaces, characterized by fixed values of ~ and C and by 
increasing Weber numbers, are represented in fig. 7, the inner surfaces corresponding to the 
lower Weber numbers. As it will be seen, at a certain distance from the atomizer outlet, the jet 
surface becomes cylindrical, the centrifugal forces due to the v, component being balanced by 
the surface tension. Actually, either the jet breaks into droplets before becoming cylindrical 
(and the droplets successively move on a conial surface) or the drag of the gas into which the 
jet is discharged (in the case of small v, components) causes it to converge towards the axis, 
forming a bubble. 


In order to control the accuracy of the theoretical results in comparison with the experi- 
mental ones, this theory has been applied to calculate the cone angle of the jet issuing from 
an atomizer tested by Soehngen and Grigull (atomizer 5 Ii). The calculation of and C has 
been carried out on the basis of the experiments of Doumas and Laster}. 


1 L. Doumas, Chemical.Engineering Progress, 49, 10 (1953) p. 518. 
24 
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In fig. 8 a comparison is made between the calculated and the experimental values. The 
agreement is good, the theoretical values, not taking into account the drag of the air entrained 
by the jet, being greater than the experimental ones. The difference is increasing for the higher 
rates of flow, where the drag is likely to be more pronounced. 

In this particular comparison, the Weber numbers being rather high, the liquid jet is unbroken 
only for a short stretch from the atomizer outlet, and the measured cone angle refers to the 


| theory 


‘experiments 


Fi V [em’sec| 
0 a 4 6 8 10 1 14 16 8 20 


Fig. 8. Comparison between experimental and calculated values of the liquid cone angle. 


30 


surface made by droplets. No difficulties have however been met in assessing the theoretical 
angle ¢, because the calculated surface is very nearly a cone at a short distence from the atomizer 
outlet. For lower Weber numbers, however, in order to compare theoretical with experimental 
values, itis necessary to evaluate the distance, from the atomizer, where the liquid surface breaks 
into droplets. Experimental values on this distance, in dimensionless form, are to be found 


in a paper by Weinberg}. 


(Eingegangen am 11. Februar 1956.) 
Anschrift des Verfassers: Prof. L. J. Valdenazzi, Genova (Italien) University of Engineering Via Montallegro 1 


* Weinberg, Proc. Inst. Mech. Eng., 1B, (1952—53), p. 240. 
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Stationdre Schwingungen in nichtlinearen dynamischen Systemen 
mit Totzeiten * 


Von K. Magnus 


1. Einfithrung. In Regelkreisen treten haufig Glieder auf, die mit Totzeiten behaftet sind. 
Bei ihnen verstreicht aus irgendwelchen geratetechnischen Griinden stets eine gewisse Zeit 
— die Totzeit, auch Laufzeit genannt —, bis sich eine am Eingang des Gliedes auftretende 
Ursache als Wirkung im Ausgang bemerkbar macht. Soll beispielsweise das Mischungsverhaltnis 
zweier Fliissigkeiten oder Gase auf Grund von Messungen an der Mischung geregelt werden, so 
darf die Messung selbst erst nach einer gewissen Mindestmischstrecke erfolgen. Die zum Durch- 
laufen der Mischstrecke gebrauchte Zeit geht als Totzeit in den Regelvorgang ein. Bei Ver- 
wendung von elektrischen Relaiseinrichtungen verstreicht zwischen dem Einschalten des Steuer- 
stromes und der Kontaktgabe eine Ansprechzeit, die als Totzeit zu betrachten ist. 

Mathematisch macht sich das Auftreten von Totzeiten dadurch bemerkbar, da® in den 
Differentialgleichungen Glieder vorkommen, bei denen das Argument um den Betrag der Totzeit 
in der Zeitachse verschoben ist. Die Gleichungen bilden dann ein System von einfachen Funk- 
tional-Differentialgleichungen, deren Lésungsmannigfaltigkeit im allgemeinen viel gréBer ist 
als bei demselben System ohne Totzeiten. Fiir lineare Systeme mit verzégerten Argumenten 
liegen zahlreiche Untersuchungen vor. Perron! und E. Schmidt? untersuchten Systeme n-ter 
Ordnung und gewinnen durch Diskussion der Wurzeln einer transzendenten charakteristischen 
Gleichung Bedingungen fiir den allgemeinen Charakter der vorhandenen Lésungen. Collatz* 
behandelt ein System zweiter Ordnung und gibt Stabilitaétsbereiche in Abhangigkeit von den 
Systemparametern an. Hahn? gelingt es, fiir dasselbe System die Lésungen in Form konver- 
genter Reihen abzuleiten. Die Anwendung der eleganten Verfahren der Laplace-Transformation 
zur Lésung von Totzeitproblemen findet man bei verschiedenen Autoren z. B. Bulgakow®, Kiipf- 
miiller® und Oldenbourg-Sartorius’. 

Der einzige, dem Verfasser bekannt gewordene Versuch, auch nichtlineare Systeme mit 
Totzeitgliedern zu behandeln, findet sich in einer Arbeit von Popow®, die sich mit der Anwendung 
des von Krylow und Bogoljubow angegebenen Verfahrens der ,,Harmonischen Balance“ auf 
Systeme beliebiger Ordnung beschaftigt. Unter Verwendung der Nyquistschen Stabilitatskriterien 
gewinnt Popow Bestimmungsgleichungen fiir Amplitude und Frequenz von eventuell vor- 
handenen Dauerschwingungen. In einer friiheren Arbeit? konnte gezeigt werden, daf} sich durch 
Kombination des Ansatzes von Krylow-Bogoljubow mit den Routh-Hurwitzschen Stabilitats- 
kriterien ein sehr allgemein anwendbares Verfahren zur Bestimmung der wichtigsten System- 
eigenschaften ableiten laBt. Das Ziel der jetzigen Untersuchungen ist der Nachweis, dai} durch 
eine geringfiigige Modifikation des Popowschen Ansatzes nun auch Totzeitprobleme in das friiher 
gegebene Schema eingeordnet werden kénnen. Es laf t sich also auch bei Auftreten von Totzeiten 
ein lineares Ersatzsystem ableiten, das keine Funktional-Differentialgleichungen enthalt und 
zu dessen Lésung bekannte Methoden verwendet werden konnen. 

Der Ansatz der Harmonischen Balance liefert eine iiber die Méglichkeiten der Methode der 
,.Kleinen Schwingungen“ erheblich hinausgehende Theorie der ersten Naherung, die insbesondere 
dann vorteilhaft verwendet werden kann, wenn in den zu untersuchenden Systemen stationadre 


* Diese Untersuchungen wurden im Rahmen eines Forschungsprogrammes durchgefiihrt, fiir dessen 
Unterstiitzung der Deutschen Forschungsgemeinschaft zu danken ist. 
1 Q. Perron, Math. Z. 45 (1939), S. 127. 
E. Schmidt, Math. Ann. 70 (1911), S. 499. 
L. Collatz, Z. angew. Math. Mech. 25/27 (1947), S. 60. 
W. Hahn, Z. angew. Math. Mech. 34 (1954), S. 316. 
B. W. Bulgakow, Schwingungen, Bd. I, Kap. 4. Moskau 1949. 
K. Kiipfmiiller, Tagungsheft ,,Die Laplace-Transformation und ihre Anwendung in der Regelungs- 
technik**. Munchen 1955. 
7 R.C. Oldenbourg und H. Sartorius, Dynamik selbsttatiger Regelungen. Miinchen 1944. 
8 E. P. Popow, Doklady Akad. Nauk SSSR. 98 (1954), S. 545. 
9 K. Magnus, VDI-Forschungsheft Nr. 451. Diisseldorf 1955. 
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Schwingungen auftreten kénnen. Die im allgemeinen nichtharmonischen Schwingungen werden 
bei diesem Verfahren durch die Grundharmonische ihrer Fourierzerlegung ersetzt. Die Fehler 
des Verfahrens werden um so geringer scin, je weniger die Schwingungen von der Sinusform 
abweichen. Allgemeingiiltige Fehlerabschatzungen sind zwar bisher noch nicht gelungen, doch 
zeigte die Berechnung einiger auch streng zuganglicher Sonderfalle, da die Fehler meist er- 
staunlich gering sind, selbst dann, wenn wegen starker Abweichungen der Schwingungen von 
der Sinusform die Brauchbarkeit des Naherungsansatzes von vornherein nicht feststeht. Wenn 
in einem System harmonische Schwingungen méglich sind — wie es im allgemeinen bei linearen 
Totzeitproblemen der Fall ist —, dann ist auch die Lésung des linearen nicht funktionalen Ersatz- 
systems streng giiltig. 

2. Herstellung eines linearen Ersatzsystems. Wie in der bereits erwahnten friiheren Arbeit 
gezeigt wurde, kann man fiir ein System ohne Totzeit, das durch die Differentialgleichungen 


4, = Sf) (i =21, 25 eee 


mit den nichtlinearen Funktionen f;, beschrieben wird, ein lineares Ersatzsystem 
n 
Pea >! (Gi) 0) ar Ree cs 70} (2) 
ot 


angeben, dessen Koeffizienten durch eine Integraltransformation aus den nichtlinearen Funktionen 
gewonnen werden: 


2x 
1 : : 
ci = ig | Sins sing) sin dp = K,{fi,} » (3) 
i 
2a 
* 1 . 
oe ara | fin Ay sing) cos yp dp = K, {fi,}. (4) 


Diese kurz als K-Transformation bezeichnete Operation verwandelt die Funktionen f;, der Vari- 
ablen x, in die Funktionen a;, bzw. aj, der Variablen A,. Die Transformationsformeln (3) und (4) 
werden durch Einsetzen des harmonischen Ansatzes 
x, = A, sin (wt — y,) (5) 
in die beiden Systeme (1) und (2) gewonnen. Bei harmonischer Veranderung der Argumente x, 
werden auch die auf den rechten Seiten von (1) stehenden Funktionen f;,(x,) periodische Funk- 
tionen der Zeit sein, die in Fourierreihen zerlegt werden kénnen. Bericksichtigt man von dieser 
Zerlegung nur jeweils das erste Glied, also die Grundharmonische, so kommt man durch einen 
Koeffizientenvergleich der rechten Seiten von (1) und (2) zu den angegebenen Transformations- 
formeln (3) und (4). 
Dieser Ansatz kann nun in einfacher Weise auch fiir Systeme mit verzégerten Argumenten 


erweitert werden. Wenn die einzelnen Argumente x, die Totzeiten t, haben, so hat man an 
Stelle von (1) das System: 


Se > Pulpit ¢ 2D Sone 


Als Ersatzsystem wird auch jetzt wieder der frithere Ansatz (2) gewahlt, nur sind fiir die Koeffi- 
zienten an Stelle von (3) und (4) andere Bestimmungsgleichungen abzuleiten. Geht man mit 
dem Ansatz (5) in die rechten Seiten von (6) ein, so erhalt man wieder periodische Funktionen 
der Zeit, die sich in eine Fourierreihe zerlegen lassen. Bei Vernachlassigung der héheren Har- 
monischen dieser Entwicklungen ergibt der Koeffizientenvergleich mit den rechten Seiten von (2) 
die Bestimmungsgleichungen 


1 
iy = | fi, [A, sin w (t — T,)] sin wt d{wt) , (7) 
0 


2x 
1 ; 
Cn SER lic [A, sin w(t —T,)] cos wt d(at) . (8) 
0 
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Setzt man darin 
sin wt = sin [w (t—T,) + @ T,] 
= sin w (t — T,) cos w t, + cos w (t —T,) sing T,, 


so folgt unmittelbar aus (7) 


2x 


1 i ; 
i a. fi, [A, sin w (t—T,)] sin @ (t — r,) cos wt, d(wt) 
0 


2x 


1 : : : 
+ aA [a , [A, sin w (t—T,)] sin w (t—t,) sinw 1, d(wt) 
i 


oder 
a;, = cos wt, K, {fi,} + w sin t, K, {fi,} ; o) 


denn bei der Ausfiihrung der Integration ist die Verschiebung des Integrationsbereiches um den 
Betrag wt, ohne Belang. Entsprechend folgt durch Einsetzen von 


- cos Wt = cos w(t —T,) cos wT, — sin w(t —T,) sin wt, 


in (8) 
2x 
1 : 
at, = oe Ar | fi» [A, sin @ (t —Tt,)] cos w(t —Tt,) cos w 7, d(wt) 
0 
aa Oe ee ; 
-|- Fe wes fi, [A, sin @ (t — t,) sin w (t —T,) sin w 7, d(ot) 
0 
oder 
at = cos w 1, Ke { fiy} — ——~" K, { fiy} - (10) 


Die Koeffizienten des linearen Ersatzsystems (2) sind nunmehr grundsatzlich frequenzabhingig, 
wahrend eine Frequenzabhangigkeit bei Systemen ohne Totzeiten nur dann auftrat, wenn Funk- 
tionen f;, vom Hysteresetyp vorhanden waren. 

Das lineare Ersatzsystem ergibt in bekannter Weise eine charakteristische Gleichung 


(Opp p> —0;, A | Cue” 0, (11) 
v= 0 
deren Koeffizienten c, fiir die weitere Rechnung gebraucht werden. 6;,ist das bekannte Kronecker- 
symbol (6;;= 1, 6;, = 0 fiir 1 + 7). 

3. Dauerschwingungen und ihre Stabilitét. Dauerschwingungen treten auf, wenn die Koeffi- 
zienten von (2) solche Werte annehmen, da sich das System gerade auf der Stabilitatsgrenze 
befindet. Legt man die Routh-Hurwitzschen Stabilitatskriterien zugrunde, so muf8 H,>0 
i —1,-% ..., n—2), aber 

tie — 0) (12) 
sein!, R ist die Routhsche Diskriminante, H, sind die Hurwitzdeterminanten. Die auf der 
Stabilitatsgrenze geltenden Frequenzen w, bekommt man in einfacher Weise aus 

oi = tes (Hy =1; nS 3), (13) 

n—2 
wobei c, das konstante Glied der charakteristischen Gleichung (11) ist?. Die Beziehungen (12) 
und (13) enthalten als Unbekannte die n+1 Groen m, und A,. Will man diese bestimmen, 
so sind also weitere n—1 Gleichungen notwendig, die sich aus dem linearen Ersatzsystem 


1 Vgl. FuBnote 9 von S. 341, dort Kap. 2.2. 
2 L. Cremer, Z. angew. Math. Mech. 25/27 (1947), S. 161. 
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unschwer ableiten lassen. Es sind die n—1 voneinander unabhangigen Verhiltnisse der Ampli- 
tuden A,, die wie folgt gewonnen werden kénnen: 


Mit dem Ansatz 


folgt aus (2) 


AWS (ayedeeeae ae (14) 


v=1 


Wiahlt man A, als Bezugsamplitude und bildet die Amplitudenverhaltnisse 


geen’ 


9 
A, 


so folgt ein inhomogenes lineares Gleichungssystem zur Bestimmung der n— 1 GroBen %;: 


n=1 


DS x, (a, + af, 4 — 0;,a) =— ain — aid (i =T, 2, aa) n—1l). (15) 
vy=1 } 
Fiir 2 ist dabei pepe 
A=] Ws (j=y—1), 
mit m, aus (13) einzusetzen. 
Die Lésung von (15) ist 
D 
x, =— = Mere == IIe 16) | 
Pdi: (yis-1, 2; , n—1) (16) | 


mit der Nennerdeterminante 
. * * 
Dy aaa |a;, SF OD Chin 67,)| . 


Die Zahlerdeterminante D, folgt aus Dy, wenn dort an Stelle der y. Spalte der Wert (— din —J4jn Ws) 
eingesetzt wird. Die Indizes t und y laufen jeweils von 1 bis n — 1. 

Die Werte x, fiir die Amplitudenverhaltnisse sind komplex, d. h. zwischen den Maximal- 
amplituden der einzelnen Koordinaten existieren irgendwelche. Phasenverschiebungen, die aus 
dem Verhaltnis von Imaginarteil zum Realteil leicht ausgerechnet werden kénnten. Sie sind in 
unserem Falle jedoch unwichtig, da die Phasenwinkel bei der Durchfiihrung der K-Trans- 
formationen (3) und (4) lediglich eine Verschiebung des Integrationsbereiches verursachen, also 
auf das Ergebnis keinen EinfluB haben. Wir interessieren uns daher weiterhin nur fiir die Betrage 
der Amplitudenverhaltnisse und setzen 


|4,| == |] 4. (v= 1, 2,...,n—1), (17) 


wobei A, als Bezugsamplitude natiirlich reell angenommen werden kann. 

Mit (12), (13) und (17) hat man nunmehr n+1 Gleichungen zur Bestimmung der n+1 
Unbekannten, also der Frequenz und der Amplitudenverteilung von Dauerschwingungen zur 
Verfiigung. 

Um auch die Stabilitat der Dauerschwingungen beurteilen oder die Frage beantworten zu 
kénnen, ob ein gegebenes System bzw. seine Stabilitatsgrenze gefahrlich oder ungefahrlich ist 
(im Sinne des von Bautin eingefiihrten Gefahrlichkeitsbegriffes), mu® auch das Verhalten in der 
Nachbarschaft der Stabilitatsgrenzen untersucht werden. In dieser Nachbarschaft ist (12) nicht 
mehr erfillt. Auch die Frequenz (13) gilt nicht mehr streng, jedoch wird dieser Wert fiir die 
weitere Rechnung beibehalten. Es wird also fiir Amplitudenwerte, die nicht den oben aus- 
gerechneten Werten von Dauerschwingungen entsprechen, fiir die Frequenz derjenige Wert ein- 
gesetzt, der sich ergeben wiirde, wenn die gewahlte Amplitude zu einer Dauerschwingung ge- 
hérte, wenn sich also das System genau auf der Stabilitatsgrenze befinden wurde. Dieser Nahe- 
rungsansatz macht die weitere Rechnung erst méglich. Er hat zur Folge, daB die fiir eine Um- 
gebung der Stabilitatsgrenze erhaltenen Aussagen quantitativ etwas verfalscht werden, jedoch 
werden Fragen qualitativen Charakters wie z.B. die nach der Stabilitat und der Gefahrlichkeit 
in richtiger Weise beantwortet. 

Um das Verhalten eines Systems in der Umgebung der Stabilititsgrenze untersuchen zu 
kénnen, miissen Veranderungen der Amplituden A, betrachtet werden. Da die Koeffizienten a;, 
und a}, nach (3) und (4) Funktionen der A, sind, so kann man (17) auch als n—1 Beziehungen 
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zwischen den Amplituden A, auffassen: 


AA a0 6) 1, 2...., n—1), (18) 


wobei (13) zu beriicksichtigen ist, da auch @, in (17) vorkommt. 
Die Lésungen des Gleichungssystems (18) seien 


ay = §$) (A,) (» =I, 2,..., n—1) r (19) 


Zu jedem Wert der Bezugsamplitude A, gehért also eine Schar von Werten A, bis 4,_1. Folglich 
kénnen auch zu jedem A, die Werte der Koeffizienten a;, und a*, angegeben werden. Bei einer 
Darstellung im Raume der Koeffizienten gibt es also zu jedem A, einen Punkt, der sich bei 
Veranderungen von A, langs einer Kurve, der ,,4-Kurve‘‘ bewegt. In friiheren Untersuchungen 
wurde nun gezeigt, daB aus der Lage dieser A-Kurve zur Stabilitatsgrenzkurve oder Stabilitats- 
grenzflache die interessierenden Eigenschaften von Systemen, deren Ersatzsystem die Form (2) 
hat, abgelesen werden kénnen. Die Ergebnisse kénnen auch analytisch formuliert werden, wenn 
man die Abhangigkeit der die Stabilitatsgrenze bestimmenden Routhschen Diskriminante von 
der Bezugsamplitude A, =A beriicksichtigt. R ist nach Definition eine Funktion der Koeffi- 
zienten a;, und aj,. Da diese von den A, abhangen und wegen (19) ist also R auch eine Funktion 
von A, deren Differentialquotient 


, dR 1 wo /aR da, | dR %*,\ a4 
= _—a=_————__ = } Z 
P= Seay Se ee 0) 


alle wichtigen Eigenschaften zu bestimmen gestattet. Es gilt!: 


1. Die aus (12), (13) und (17) ausgerechneten Dauerschwingungen sind stabil, wenn 


(R’)rao>0 . (21) 
ist. Ist dieser Wert negativ, so sind die Dauerschwingungen instabil. 


2. Eine Stabilitatsgrenze ist gefahrlich im Kleinen (im Sinne von Bautin), wenn in ihrer Nahe 
(R40 <0 (22) 
gilt. Bei Ungefahrlichkeit im Kleinen ist der gleiche Ausdruck positiv. 
3. Eine Stabilitatsgrenze ist gefahrlich im GroBen, wenn in ihrer Nahe 


(R))4-n <0 (23) 
gilt. Bei Ungefahrlichkeit im Grofen ist der Ausdruck positiv. 


4. Die Erregung von Dauerschwingungen ist weich, wenn R(0)<0 ist und zwischen dem 
Amplitudenwert der betrachteten Dauerschwingung A = Ay und dem Wert 4A =0 keine weitere 
Lésung des Systems (12), (13), (17), also keine weitere mégliche Dauerschwingung existiert. 
Hieraus folgt z. B. sofort, daB Systeme, die nur eine, und zwar stabile Dauerschwingungsamplitude 
besitzen, stets weich erregt werden kénnen. 


5. Die Erregung von Dauerschwingungen ist hart, wenn auBer der betrachteten Amplitude 
A = Aj der Dauerschwingung noch mindestens ein kleinerer Wert A = A, < Ay existiert, der zu 
einer instabilen Dauerschwingung gehért. 

Beziiglich weiterer Einzelheiten mu auf die zitierte friihere Arbeit? verwiesen werden. 


4. Beispiele. Es sollen im folgenden zwei besonders einfache Beispiele behandelt werden, die 
auch streng mit elementaren Mitteln losbar sind. Der Fehler der Naherungsrechnung ergibt sich 
dann unmittelbar durch Vergleich. Ein komplizierteres nicht mehr elementar lésbares Problem 
dieser Art soll gesondert veréffentlicht werden. 

a) Ein System erster Ordnung. E. Schmidt® hat nachgewiesen, daf} die lineare Funk- 
tional-Differentialgleichung 


x +ex(t—t) =0, (24) 
harmonische Liésungen von der Form 
x=acoswt+ bsinwt, . (25) 
1 Vel. FuBnote 9 von S. 341, dort Kap. 2.3. 


2 Vel. FuBnote 9 von S. 341. 
3 Vgl. FuBnote 2 von S. 341. 
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mit beliebigen Konstanten a und b unter der Bedingung besitzt, daf 


@ =o nee (26) 


Tt 


ist. In jedem anderen Falle gibt es nur die triviale Lésung x = 0. | 
Nach dem hier beschriebenen Verfahren bekommt man fiir das Totzeitglied einen linearen 


Ersatzausdruck 
ex(¢—t) =ax+a*x. (27) 
Mit f(x) = cx folgt aus (3) und (4) 
K, {fj} =e, KAS j= 0% 


so da& man aus (9) und (10) die Ersatzkoeffizienten 


: csin WT 
d= € C08 ®t, Ct) eae (28) 


erhalt, Fir (24) ergibt sich somit die lineare gewohnliche Ersatz-Differentialgleichung 


#(1 922") + xe cos7 =0. (29) | 


Wenn periodische Lésungen dieser Gleichung mdglich sein sollen, muB notwendigerweise 


csin ® T 


ieee =0, ccosat=0O0, (30) | 
gelten. Daraus folgt aber die schon von Schmidt streng abgeleitete Beziehung (26). Die Fre- 
quenz wird also von der Naherungsliésung fehlerfrei wiedergegeben. Die Amplitude bleibt un- 
bestimmt, wie es bei einem linearen Problem nicht anders zu erwarten ist. 

Nun soll (24) durch Einfiihrung der Vorzeichenfunktion sgn zu einer nichtlinearen Differen- 
tialgleichung gemacht werden: 


x + esgn x(t—t) =0. (31) | 
Es liegt hier ein besonders einfacher Fall einer ,,Zweipunktregelung‘‘ vor, da das zweite Glied 
nur der beiden diskreten Werte -+-c und —c fahig ist. Die strenge Lésung von (31) ist 


x=+e, x= tet+ %. 
Beginnt man die Zeitzéhlung beim Nulldurch- 
gang der x-Koordinate, so bekommt man als 


x(t)-Kurve eine Dreieckschwingung, wie sie 
Abb.1 zeigt. Ihre Ampltiude ist 


Al = OB 5 (32) 
ihre Frequenz 
I 


; Nach der Naherungsrechn i 
Abb. 1. Amplituden-Zeit-Kurve fiir die periodische Grundlésung ) cAnUne hat man zunachst 


von (31), (¢ = Totzeit), esgn x(t—t) =ax+a*x. (34) | 


Fir die Sprungfunktion f(«) =c sgn x folgt aus (3) und (4) bzw. aus einer Tafel von K-Trans- 
formierten: 


K{f}=25, KAf}=0. 


1 


Dann ergibt (9) und (10) 


a ns ee a* chery 
= = 5 SS i 
aA ; mAw a Os 


so daB als lineare Ersatzdifferentialgleichung fiir (31) 


: Ave! tar Ac 
x (1 Sa e r) “+ x4 os wt =0 (35) 
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erhalten wird. Periedisehe Lésungen dieser Differentialgleichung sind nur méglich fiir 


ATC 4c 
DS A Oe ag COs ot= 0. (36) 


Da die Amplitude 4 nur positiver Werte fahig ist, folgen also Lésungen mit den Frequenzen und 
Amplituden 


o =1(%4 20a] ah, eae (37) 


4c Act 
Gay a . 
n(F + 2na] 


Im Faille n =0 ist der Frequenzausdruck identisch mit der strengen Lésung (33), wahrend bei 
der Amplitude 4 der Faktor 8/7? an Stelle von 1 steht, also ein um 18,9° zu kleiner Wert erhalten 
wird. 

Fir andere Werte von n bekommt man Schwingungen, bei denen innerhalb der Totzeit ein 
mehrmaliges Umschalten erfolgt. Wie man sich leicht iiberlegt, folgen dafiir aus (37) jeweils 
die streng richtigen Frequenzwerte, jedoch Amplitudenwerte, die stets um den Faktor 8/z? zu 
klein sind. Derartige héherfrequente Schwingungen sind zwar theoretisch und unter gewissen 
Voraussetzungen auch praktisch méglich, jedoch geniigt zu ihrem Zustandekommen nicht die 
Vorgabe irgendwelcher bestimmter Anfangsbedingungen; es muf} vielmehr eine Vorgeschichte 
von der Zeitdauer tT mit genau definierten Nulldurchgingen vorhanden sein. 


b) Eine unstetige Regelung zweiter Ordnung. Erganzt man die Differential- 
gleichung (31) noch durch ein ,,Schleppglied‘‘, so findet man 


Tx +x +csgnx(t—rt)=0, (28) 


wobei T als eine Zeitkonstante aufgefaBt werden kann. Eine Differentialgleichung dieser Art 
erhalt man zum Beispiel bei der naherungsweisen Berechnung der Kursregelung von Schiffen 
oder Flugzeugen, wenn nur die Momentengleichung beriicksichtigt wird und der Regler selbst 
nach dem Schwarz- Wei-Prinzip (Zweipunktregler) arbeitet. 

Zunachst soll eine periodische strenge Lésung von (38) angegeben werden. Im Bereich 
x(t—t) >0 folgt aus (38) 


€ 
Bier ey | 
(39) 


t 


e=—BTe Tct+B,| 


mit den beiden Integrationskonstanten B, und B,. Fur den Bereich x (t—t) < 0 erhalt man 
entsprechende Ausdriicke, bei denen nur das Vorzeichen von c verdndert ist. Wenn sich die 
Lésungen aus beiden Bereichen zu einer periodischen Gesamtlésung zusammensetzen sollen, so 
miissen die folgenden Anstiickelungsbedingungen gelten: 


x(t;) = — x(0), 
x(t,) = — x(0) , (40) 
x(t, —t) =0. 


Die ersten beiden Bedingungen bringen zum Ausdruck, daf periodische Lésungen in beiden 
Bereichen symmetrisch sein miissen, di¢ dritte Bedingung definiert die Zeit t;, in der die Lésung (39) 
die Grenze ihres Giiltigkeitsbereiches erreicht. Hierbei ist t, gleich der halben Schwingungszeit. 

Durch Einsetzen von (39) in (40) bekommt man zuniachst zwei Gleichungen fiir B, und B,, 
aus denen die beiden Integrationskonstanten ausgerechnet werden kénnen, 


Lee . (41) 
B=tT +t. 
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sowie eine Gleichung fiir die unbekannte Schaltzeit t, 


Reels y pane pyc (42) 


a 


eee ee Ve 


Hat man aus dieser Gleichung t; — zum Beispiel durch Iteration — gefunden, so ist die Frequenz 
der Schwingung 


ie (43) 
hy 


Die Amplitude der periodischen Lésung folgt aus x(t) =0 und x(t) =A. Aus (39) findet man 
leicht 
c 
eaceisrineas (44) 


1 
A= B,—eT(1—In $), (45) 
1 


oder unter Beriicksichtigung von (41) nach Umformung: 
t, 
Ae . + Tn Fe 4 “7 )h (46) 


Fiir gegebene Parameter T, t und c des Systems bestimmt man also zunachst aus (42) den 
Wert t,, der wegen (43) unmittelbar die Frequenz der Schwingung ergibt. Danach kann aus (46) 
der zugehérige Wert fiir die Amplitude A errechnet werden. Man beachte, da die Frequenz 
der Schwingung nicht von dem Parameter c abhangt, wahrend die Amplitude proportional zu 
ihm ist. 

Nach der Naherungsrechnung bekommt man zunichst denselben linearen Ersatzausdruck (34) 
wie im vorigen Beispiel. Auch die Ausdriicke fiir die Koeffizienten a und a* andern sich nicht. 
Die lineare Ersatzdifferentialgleichung wird demnach 


c 


s 4 
8 Ale 


Z ; de 

sin WT) *X — cos Oe AT 
peer: rn 
Die Routhsche Diskriminante reduziert sich in diesem einfachen Fall auf den Faktor von x? 
wahrend das Quadrat der Frequenz in bekannter Weise als Quotient der Faktoren von x und x 
erhalten wird. Man hat somit zur Bestimmung stationarer Schwingungen die beiden Gleichungen 


[ke == Il = te snwot=—QO0, (48) 
azAw 

9. . HE p A 

cage Gree as ee 


Durch Elimination der Amplitude folgt daraus eine Gleichung fiir die Frequenz: 
M 
tgot= TF: (50) 


Da A nur positiver Werte fahig ist, so sieht man aus (48) und (49), daB® nur solche Lésungen 
fir w physikalisch sinnvoll sind, fiir die coswm t > 0 und sinwt > 0 ist. 
Hat man w ermittelt, so folgt die Amplitude sofort z.B. aus (48) 


A Cuae 
A = Fp not. (51) 


Wie bei der strengen Lésung ist A proportional zu c, wihrend @ von c unabhiangig ist. 

Fir ein festes Wertepaar c und T ist der prozentuale Fehler der Naherungslésungen fur @ 
und A als Funktion der Totzeit t ausgerechnet und in Abb. 2 aufgetragen worden. 

Zur Untersuchung der Stabilitat sowie der Gefahrlichkeit wird gebildet 


dR aR , @Rdo 


dA gue eon (52) 
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Mit dem Ausdruck von (48) bekommt man 


dR Aker ees w T (sin t —T cos wT) 
= = —— sin WT — 
dA nA? a A (2 coswt + tSin wT)" 


(53) 


Nach (21) entscheidet das Vorzeichen dieses Ausdruckes fiir den Fall R = 0 iiber die Stabili- 
tat. Beriicksichtigt man (50) und (51), so folgt aus (53) 


a _ T+rt+te?T 
dAjp=6) AQT Ez) 
Die ausgerechnete Dauerschwingung ist also stabil. 


Zur Kladrung der Frage der Gefahrlichkeit im Kleinen ist in (53) der Grenziibergang A —> 0 
vorzunehmen. Wegen cos wt—> 0 und sinwt—> + 1 findet man 


= 0 (54) 


(a) 
we + pes ae CO) < 
Die Stabilitatsgrenze des Systems ist also nach (22) ungefahrlich im Kleinen. 


Im anderen Grenzfalle A — oo 
bekommt man 


(34 0 
Tale a oue 


%-Fehler 


10 
A 
5 
@ 
Ci 
0 O5 7sek 


Abb. 2. Prozentualer Fehler der Naherungen fir die 
Amplitude A und die Frequenz » als Funktionen 
der Totzeit 1 fiir das Beispiel b). Abb. 3. Stabilitatsdiagramm fir das Beispiel b). 


webei, wie man leicht feststellt, der Differentialquotient von positiven Werten gegen Null geht. 
Obwohl die Stabilitatsgrenze des Systems nach dem Kriterium (23) weder als gefahrlich noch als 
ungefahrlich im GroBen bezeichnet werden kann, wird sie praktisch ungefahrlich im Grofen sein, 
da bei groBen, aber endlichen Werten fiir A auf jeden Fall 
dR 
ec 
ist, : 

Man kann die hier analytisch aus den Kriterien (21), (22) und (23) ermittelten Eigenschaften 
auch durch eine Untersuchung des Verlaufs der ,,A-Kurve“ bekommen. Wir zeichnen diese Kurve 
zusammen mit der Stabilitatsgrenzkurve in der Ebene der beiden Koeffizienten b und 6*, die 
- durch Umschreiben von (47) in die Form 


x+t+bx + b*x=0 


erhalten werden. 6 ist dann bis auf den Faktor 1/T gleich dem Wert R von (48), b* ist gleich 
dem Wert w? von (49). Aus (49) kann nun zu jedem Wert von A der zugehérige Wert von @ aus- 
gerechnet werden. Dann folgt aus (48) fiir jedes Wertepaar A, « der entsprechende Wert von R 
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und damit auch von b. Die so erhaltene A-Kurve in der 6b*-Ebene zeigt Abb. 3, wobei 
als Parameter auf der Kurve die Werte von A angeschrieben wurden. Die Kurve kommt aus 
dem negativ Unendlichen und geht mit A — oo gegen den Punkt b = 1; b* = 0. Die Stabilitats- 
grenzkurve (R-Kurve) ist nach (48) durch R = 0 gegeben, d. h. in unserem Falle durch b = 0. 
Im oberen Quadranten ist R > 0 (stabiler Bereich), im unteren ist R < 0 (instabiler Bereich). 

A-Kurve und R-Kurve schneiden sich im Punkte A, = 1,12; b* = w? = 0,740. Damit 
sind Amplitude und Frequenz der Dauerschwingung bestimmt. Ihre Stabilitat folgt einfach aus 
der Tatsache, daB die A-Kurve mit wachsendem A vom instabilen Bereich kommend in den 
stabilen eindringt. Fiir kleine Werte von A lauft die A-Kurve auf den Stabilitatsbereich zu, 
folglich ist die Stabilitatsgrenze als ungefaihrlich im Kleinen anzusprechen. Im Grenzfalle 
A-> o wird die A-Kurve parallel zur Stabilitatsgrenze; das System ist also streng genommen 
im Gro®en neutral, jedoch kann es praktisch als im Groen ungefahrlich angesehen werden, da 
die A-Kurve fir endliche groBe Werte von A noch immer von der Stabilitatsgrenze fortlauft. 
Nach den im Kap. 3 angegebenen Kriterien kénnen die Dauerschwingungen weich erregt werden, 
denn es existiert nur ein Schnittpunkt von A- und R-Kurve, und dieser geh6ért zu einer stabilen 
Dauerschwingung. 


(Eingegangen am 14. Februar 1956.) 
Anschrift des Verfassers: Dr. Kurt Magnus, Freiburg/Breisgau, Unterlinden 7. 


Konstruktive Grundziige und praktische Erfahrungen 
beim Bau und Betrieb von Stahlwasserbauten 


Von Dipl.-Ing. Franz Kohler, Miinchen. Mit 75 Abbildungen. VII, 135 Seiten Gr.-8°. 1956. 
Steif geheftet DM 13,20 


Inhaltsiibersicht: Erster Teil: Wehrverschliisse. A. Allgemeines: Bewegliche Wehrverschliisse. An- 
triebe. Gewichte und Kosten. — -B. Dichtungen: Allgemeine Anforderungen. Ausfiihrungsarten. — 


C. Wehrbeheizung. — D, Unterwasserrostschutz. — E. Berechnung von Stahlwasserbauten: Stahlkon- 


struktionen. Antriebe. — Zweiter Teil: Schleusenversehliisse. A. Allgemeine Aufgaben der Fiill- und 
Ee -Entleerungsorgane: Binnenschleusen. Seeschleusen. — B. Statische und konstruktive Grundsitze: All- 
gemeinés. Lastannahme. Statische, konstruktive und bauliche Besonderheiten. — C, Antriebe: All- 


' gemeines. Huborgane. Motorleistungen. Bemessungsgrundlagen. Anwendung bei verschiedenen Ver- 


schluBsysteméen. — D. Unterwasserrostschutz. — E. Zusammenfassung. — Sachverzeichnis. 
4 


Das Fachgebiet des Stahlwasserbaues hat durch die nach Kriegsende wieder einsetzende Bautitigkeit 
& neuen Aufschwung erfahren. Neben den im Wehr- und Schleusenbau allgemein auftretenden Problemen 
bautechnischer Art haben vor allem auch die Verschlu8konstruktionen und die Antriebe weaiteetehe Ver- 
besserungen aufzuweisen. Teilweise ist auch eine neue Entwicklung SaLAabe MaBgeblich beteiligt an 
diesen Fortschritten sind die Erfahrungen, die in den letzten 30 Jahren beim Bau und Betrieb von Stahl- 
wasserbauten gesammelt werden konnten. Diese zu vermitteln sind den gegenwartigen Stand der Ent- 


wicklung aufzuzeigen, ist der Zweck dieser zusammenfassenden Abhandlung. 


—~ 


Veréffentlichungen zur Erforschung der DruckstoB- 
probleme in Wasserkraftanlagen und Rohrleitungen 


Herausgegeben von Dr.-Ing. Friedrich Télke, Docteur és Sciences h. c., 0. Professor an der Technischen 
Hochschule Stuttgart. s 


Zweites Heft: Mit 225 Abbildungen. V, 149 Seiten 4°. 1956. _ DM 37,50 


Inhaltsiibersicht: Uber die konstruktive Gestaltung der druckstoSgefahrdeten Teile von Wasserkraft- 
anlagen. Von Professor Dr.-Ing. Fr. Télke, Stuttgart. — Experimentelle Untersuchungen tber die Stiitz- 
wirkung der Druckschachtpanzerung der Druckrohrleitung eines Alpen-Wasserkraftwerkes. Von Professor 
Dr.-Ing. W. Pelikan, Stuttgart. — Nichtstationdre Strémungen in Unterwasserstollen. Von Dr.-Ing. 
H. Blind, Miinchen. — Drucksto8 und Wasserschlag beim Abschluf von Hauswasserleitungen. Von 
-Oberbaurat W. Gandenberger, Stuttgart. — Berechnung der Grofe der Druckwindkessel bei Wasser- 


werken. Von-Professor Dr. phil. K. Ludwig + und Baurat Dr.-Ing. H. Stack, Hannover. 
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Die Statik im Stahlbetonbau 


. Ein Lehr- und Handbuch der Baustatik. Von Dr.-Ing. Kurt Beyer, chem. ord. Professor an. der 


Technischen Hochschule Dresden. Zweite, vollstandig neubearbeitete Auflage. Zweiter berichtigter — 


Neudruck, Mit 1372 Abbildungen im Text, zahlreichen Tabellen und Rechenvorschriften. XII, 804 
Seiten Gr.-8°. 1956. Ganzleinen DM 66,— 


Um der regelmaSigen Nachfrage zu geniigen, wurde ein zweiter berichtigter Neudruck hergestellt. ,,Im Gegensatz zu einer 
Reihe von Lehrbiichern iiber Statik entwickelt der Verfasser nicht eine eigene neue Methode mit dem etwaigen Nachweis 
der Anwendbarkeit auf die typischen Tragwerksformen. Vielmehr fulbt er auf dem Bestand der vorliegenden wissen- 
schaftlichen Arbeiten, aus denen er in kritischer Wiirdigung das Gemeinsame der zahlreichen verschiedenartigen Methoden 
der Baustatik herausschiilt und zu einer quasi ibergeordneten Theorie des Tragwerkes verarbeitet und einheitlich zusammen- 
faBt. Hierzu dient inshesondere die eingehende Darstellung der Grundlagen der Baustatik, die die Einleitung des Werkes 
bildet. Die ersten vier Kapitel enthalten die Theorie des Tragwerkes, und zwar: Grundlagen, das statisch bestimmte Trag- 
werk, die Forminderung des ebenen Stabzuges und die statisch unbestimmten Tragwerke. Kapitel V behandelt die An- 
wendung der Theorie auf die hauptsiachlichen Tragwerksformen, wie Durchlauftrager, Stockwerksrahmen, Rahmentrager, 
Silozellen, Bogentrager, Zweigelenk- und eingespannter Bogen, Kreisringtrager und Tragerroste. Den AbschluB bildet 
Kapitel VI mit der Behandlung der Flachentragwerke, wie Platten, Scheiben und Schalen. Mit Riicksicht auf den erheb- 
lichen Umfang des gesamten Stoffgebietes ist der Text straff gefaBt, dabei erschépfend und streng wissenschaftlich, so 
daB schon eine gewisse Mitarbeit des Lesers erforderlich ist, um den vollen Gehalt der pragnanten Ausfihrungen zu er- 
sehlieBen. Jedes Teilgebict wird zunachst theoretisch behandelt und begriindet und dann an gut gewahlten Beispielen, 
die z. T. vom theoretischen Ansatz bis zum Ergebnis voll durchgerechnet sind, erlautert und mit zahlreiehen Tabellen und 
Rechenvorschriften versehen, die die Arbeit des praktischen Ingenieurs wesentlich erleichtern und abkirzen. Durch diese 
Verbindung zwischen abstrakter Theorie und der konkreten Anwendung auf den Hinzelfall wird das Werk zum wirklichen 
Handbuch der Praxis. Fiir den wissenschaftlich arbeitenden Ingenieur geben die hinter jedem Abschnitt eingeschalteten 
Literaturhinweise die Méglichkeit zu einem yertieften Studium der Sondergebiete an. Durch die weitgehende Unterteilung 


der einzelnen Abschnitte, die straffe Ordnung und Glicderumg des Stoffes, das ausfihrliche Sachregister und Verzeichnis . 


der Reehenvorschriften und Zahlenbeispiele wird die Benutzung des Werkes sowohl zum Studium der Statik wie zur 
Lésung von Einzelfragen angenehm erleichtert. Die bewunderungswiirdige Systematik und die Einheitlichkeit im Auf- 


bau der Theorie in Verbindung mit den Anwendungsbeispielen der Praxis verbirgen dem Werk einen vag ean ee 
‘ >, Baustatik® 


und stellen es in die mpaiaenlelnt paced der deutschen Fachliteratur.** 
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Bemessungsverfahren fir Verbundtrager 


Von Dr.-Ing. Giinter Utescher, Karlsruhe. Mit 27 Abbildungen und 24 Bemessungstafeln. VII, 39 
Seiten und 24 Tafeln. Gr,-8°, 1956. : Steif geheftet DM 12,— 


/ 


Tahalteveraese tates A. Binfibning. — B. Blechtrager: Ableitung der Bemessungsgleichung. Durch- 
fiihrung der Bemessung. Der Sonderfall My = My, Kr. Der Sonderfall M,, = 0. Verstarkung des Stahl- 


.obergurtes. — C. Walzprofil-Trager: Aufstellung der Bemessungstafeln. Benutzung der Bemessungstafeln. 


Verstarkung von Walzprofilen. — D. Kriechen:-Wahl des Verfahrens. Nachweis der Kriechspannungen 
nach Fritz. Nachweis der Kriechspannungen nach Miller. — E, Sehwinden: Wahl des Verfahrens. 
Nachweis der Schwindspannungen nach Fritz. Nachweis der Schwindspannungen nach Miller. — 
F. Temperatur. G. Abstufung der Stahlirigergurte. — H. Walzprofil-Trager der Gruppe 1/DIN 4239: 
Aufstellen der Tragfahigkeitstafeln. Benutzung der Tragfahigkeitstafeln. Nachweis der Betonspannung. — 
I. Zahlenbeispiele: I. Bemessung eines Blechtrager-Verbundquerschnittes fiir Beanspruchung durch M,; 
und My. — IL. Bemessung eines Blechtrager-Verbundquerschnittes fiir Beanspruchung ‘allein durch 
My. UT. Verstarkung des Stahl-Obergurtes. — IV. Bemessung eines Walzprofil-Verbundquerschnittes 
mit Nachweis der Kriech- und Schwindspannungen. — V. Bemessung eines Verbundtragers mit Ab- 


stufung des Untergurtes. — VI. Bemessung eines Walzprofil-Verbundquerschnittes der Tragergruppe I. — 
Vereinfachter Nachweis. — Schrifttum. — Tafelanhang: Tafel 1—12 fiir Normalprofil-I-Trager, Tafel 13 
bis 24 fiir Breitflansch-I-Trager. A 


Die Zulassung der Verbundbauweise fiir den Briickenbau und Hochbau hat nicht nur neue, mit wirtschaftlichen Vorteilen , 
verbundene Konstruktionsgrundsatze gebracht, sondern auch einen erheblichen Aufwand an statischen Untersuchungen, 
Querschnittsermittlungen und Spannungsnachweisen zur Folge gehabt. Das Ziel, die zulassigen Stahl- und Beton- 
spannungen nicht nur einzuhalten. sondern auch méglichst auszunutzen, hat die Bemessung von Verbundquerschnitten 
zu einer mithsamen und zeitraubenden Probierarbeit werden lassen. / eat f 

Die vorliegende Arbeit hat sich zum Ziel gesetzt, das bisherige Ausprobieren eines mebr oder weniger willkiirlich gewahlten 
Stabltrdgerquerschnittes durch eine planmafige Querschnittsbestimmung zu ersetzen, bei der die speduupeocnanaaen 
Einfliisse des Kriechens und Schwindeng durch Naherungsformeln erfaBt werden, die in der Regel nur noch eine geringe, 
aber auch wieder planmafig Hein Fetes Korrektur erforderlich machen. } ‘ 
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